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The influence of magnetic dipole interactions (MDIs) on the magnetic properties of local-moment 
Heisenberg spin systems is investigated. A general formulation is presented for calculating the 
eigenvalues A and eigenvectors jl of the MDI tensor of the magnetic dipoles in a line (one dimension 
ID), within a circle (2D) or a sphere (3D) of radius r surrounding a given moment fli for given 
magnetic propagation vectors k for collinear and coplanar noncollinear magnetic structures on both 
Bravais and non-Bravais spin lattices. Results are calculated for collinear ordering on ID chains, 2D 
square and simple-hexagonal (triangular) Bravais lattices, 2D honeycomb and kagome non-Bravais 
lattices and 3D cubic Bravais lattices. The A and jl values are compared with previously reported 
results. Calculations for collinear ordering on 3D simple tetragonal, body-centered tetragonal, and 
stacked triangular and honeycomb lattices are presented for c/a ratios from 0.5 to 3 in both graphical 
and tabular form to facilitate comparison of experimentally determined easy axes of ordering on 
these Bravais lattices with the predictions for MDIs. Comparisons with the easy axes measured 
for several illustrative collinear antiferromagnets (AFMs) are given. The calculations are extended 
to the cycloidal noncollinear 120° AFM ordering on the triangular lattice where A is found to be 
the same as for collinear AFM ordering with the same k. The angular orientation of the ordered 
moments in the noncollinear coplanar AFM structure of GdB 4 with a distorted stacked 3D Shastry- 
Sutherland spin-lattice geometry is calculated and found to be in disagreement with experimental 
observations, indicating the presence of another source of anisotropy. Similar calculations for the 
undistorted 2D and stacked 3D Shastry-Sutherland lattices are reported. The thermodynamics of 
dipolar magnets are calculated using the Weiss molecular field theory for quantum spins, including 
the magnetic transition temperature T m and the ordered moment, magnetic heat capacity and 
anisotropic magnetic susceptibility x versus temperature T. The anisotropic Weiss temperature 0 P 
in the Curie-Weiss law for T > T m is calculated. A quantitative study of the competition between 
FM and AFM ordering on cubic Bravais lattices versus the demagnetization factor in the absence 
of FM domain effects is presented. The contributions of Heisenberg exchange interactions and of 
the MDIs to Tm and to 0 P are found to be additive, which simplifies analysis of experimental data. 
Some properties in the magnetically-ordered state versus T are presented, including the ordered 
moment and magnetic heat capacity, and for AFMs the dipolar anisotropy of the free energy and 
the perpendicular critical field. The anisotropic \ f° r dipolar AFMs is calculated both above and 
below the Neel temperature Tn and the results are illustrated for a simple tetragonal lattice with 
c/a > 1, c/a = 1 (cubic) and c/a < 1, where a change in sign of the x anisotropy is found at c/a = 1. 
Finally, following the early work of Keffer [Phys. Rev. 87, 608 (1952)], the dipolar anisotropy of 
X above Tn = 69 K of the prototype collinear Heisenberg-exchange-coupled tetragonal compound 
MnF 2 is calculated and found to be in excellent agreement with experimental single-crystal literature 
data above 130 K, where the smoothly increasing deviation of the experimental data from the theory 
on cooling from 130 K to Tn is deduced to arise from dynamic short-range collinear c-axis AFM 
ordering in this temperature range driven by the exchange interactions. 


I. INTRODUCTION 

Local magnetic moments generate magnetic dipole 
fields around them. In local-moment spin systems, the 
long-range magnetic dipole interaction between the lo¬ 
cal magnetic moments (spins) is always present. How¬ 
ever its strength is usually small compared to other in¬ 
teractions such as exchange and RKKY interactions be¬ 
tween the spins. The thermal-average magnitude of the 
interaction energy is of order E ~ p, 2 /r, where fi is 
the thermal-average value of the magnetic moment and 
r is the distance between nearest-neighbor spins. Tak¬ 
ing, e.g., )jl = 7 for Gd +3 or Eu +2 (/xb is the Bohr 
magneton) and r = 3 A gives E/k b ~ 0.02 K (kn is 
Boltzmann’s constant), which is usually very small com¬ 
pared to the other interactions between the spins. How¬ 
ever, even when the dipole interactions are weak, these 


interactions can be decisive in determining the orienta¬ 
tions of the ordered moments in magnetic structures of 
local-moment ferromagnets (FMs) or antiferromagnets 
(AFMs). 

If the distance between local moments is large enough, 
the magnetic dipole interaction can dominate the ex¬ 
change interactions in local-moment insulators and result 
in either FM or AFM dipolar ordering. Examples include 
FM ordering between Mng" 3 clusters with spin S = 12 
at the Curie temperature Tq = 0.161(2) K in mono- 
clinic Mn604Br4(Et2dbm)6, 1,2 and AFM ordering in the 
face-centered cubic (fee) diamond lattice of rare-earth R 
atoms (R = Gd, Dy, Er) in i?P 04 (Mo 03 )i 2 • 3 OH 2 O with 
Neel temperatures Tn = 0.01 K to 0.04 K. 3 

The theoretical study of magnetic dipole interactions 
and associated magnetic structures in crystals has a long 
history. In 1946 Luttinger and Tisza solved for the 
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possible magnetic structures of simple cubic (sc), body- 
centered cubic (bcc) and fee Bravais spin lattices aris¬ 
ing solely from classical magnetic dipole interactions, 
where the ordered moments all had the same magnitude 
(equal-moment magnetic structures ). 4 They found that 
the ground state for the sc lattice is an AFM state with 
propagation vector k = (|,i,0) r.l.u., whereas a FM 
state with k = 0 is the ground state for the bcc and fee 
lattices if the samples are in the shape of long thin nee¬ 
dles, but AFM structures are the most stable structures 
otherwise with k = (|, 5 ,0) r.l.u. and k = (|, i) r.l.u., 
respectively. The abbreviation r.l.u. means reciprocal lat¬ 
tice unit, where 1 r.l.u. = 2 -k/ a for cubic lattices and a 
is the cubic lattice parameter. Cohen and Keffer con¬ 
firmed using spin-wave theory that FM cannot be the 
ground state at T = 0 for pure magnetic dipole inter¬ 
actions in thin needles of a sc spin lattice but can be 
the ground state for bcc and fee lattices . 5 The magnetic 
structures of two-dimensional (2D) Bravais spin lattices 
induced by magnetic dipole interactions have also been 
investigated . 6-10 

Luttinger and Tisza also showed that in a classical cu¬ 
bic dipolar AFM in the magnetically-ordered state at 
temperature T = 0 with a magnetic field H z applied 
perpendicular to the easy axis of ordering, the compo¬ 
nent )i z of the ordered moment per spin in the direction 
of H z is proportional to H z for 0 < H z < H c and is 
equal to the saturation moment // sat for H z > H c , where 
H c is termed the critical field . 4 An expression for the 
magnetic susceptibility = fJ, z /H z for 0 < H z < H c 
was given. The high-field state with H > H c is a field- 
induced paramagnetic (PM) state in which the magnetic 
moments are ferromagnetically aligned in the direction 
of H z with )i z = n sat . According to the Weiss molecular 
field theory (MFT), precisely the same type of n z (H z ) 
behavior for the perpendicular magnetization occurs for 
both collinear and coplanar noncollincar AFMs with the 
spins interacting only by Heisenberg exchange . 11 The sus¬ 
ceptibility parallel to the easy axis at T < Tn for dipolar 
AFM ordering in a uniaxial (tetragonal or hexagonal) 
crystal has not been calculated before to our knowledge. 

The so-called pyrochlore spin lattice has attracted 
much attention over the past two decades in the con¬ 
text of spin-ice compounds . 12 This non-Bravais fee spin 
lattice with 16 spins per fee unit cell consists of a 3D net¬ 
work of corner-sharing tetrahedra formed by either the 
A or B sublattices of a pyrochlore-structure compound 
A 2 B 2 X 7 or by the B sublattice of a spinel-structure com¬ 
pound AB 2 X 4 . An example is the Ho sublattice in the 
pyrochlore compound Ho 2 Ti 2 C> 7 , where due to crystalline 
electric field effects the Ho cations behave at low T like 
Ising spins that can only point along the [ 111 ] and equiv¬ 
alent crystal directions (the Ti +4 cations are nonmag¬ 
netic). The spin-ice arrangement of the Ho moments 
at low T gives rise to a macroscopic degeneracy and a 
nonzero spin entropy at T = 0, as occurs in water ice. 
Magnetic dipole interactions between the Ho moments 
have been determined to be important to this magnetic 


behavior , 12 and hence these compounds are sometimes 
referred to as dipolar spin ices. 

On another front, dynamic magnetic fluctuations in 
long-range ordered 3D AFMs mediated by magnetic 
dipole interactions are stronger than for exchange inter¬ 
actions on the same lattice , 3,13 contrary to what might 
have been anticipated from the classical origin of the 
magnetic dipole interaction. In particular, in the cubic 
diamond lattice dipolar AFMs I?P 04 (Mo 03 )i 2 • 30H 2 O 
(R = Gd, Dy, Er), White et al. found that the suppression 
of the T —> 0 susceptibility versus 1/S 1 (S is the effective 
spin quantum number) due to quantum fluctuations was 
a factor of two stronger than predicted for the nearest- 
neighbor Heisenberg model on the diamond lattice . 3 Cor- 
ruccini and White found that within spin-wave theory, 
the 3D sc lattice exhibits quantum corrections to the 
Neel state that are also a factor of two larger than those 
of the nearest-neighbor Heisenberg AFM on the same 
lattice, indicating that dipolar magnets are more quan¬ 
tum mechanical than generally suspected, whereas the 
2D dipolar square lattice does not exhibit long-range or¬ 
der at finite temperature . 13 On the other hand, several 
authors have found that dipolar interactions in conjunc¬ 
tion with Heisenberg interactions can induce long-range 
order at finite temperatures on 2D spin lattices . 14-16 

The influence of magnetic dipole interactions on the 
magnetic properties of 3D Bravais spin lattices with unit 
cell symmetries lower than cubic has been discussed for 
particular cases. Rotter has discussed the predictions of 
dipolar interactions for the easy axis of collinear AFMs 
with AFM propagation vectors k that are determined by 
isotropic Heisenberg exchange interactions in a variety of 
collinear AFM compounds containing sc, fee, hexagonal 
and body-centered tetragonal (bet) Gd sublattices . 17 He 
found that in most cases the easy axis is consistent with 
that predicted for magnetic dipole interactions. Several 
authors calculated the local dipolar fields at a lattice site 
for general simple tetragonal and bet Bravais spin lattices 
versus a parameter not proportional to the c/a ratio . 18-20 
Maurya et al. calculated the influence of magnetic dipole 
interactions on the magnetization versus field isotherms 
of three AFMs containing Eu +2 spins-7/2 below their 
Neel temperatures of 12 K to 15 K . 21,22 

Classical Monte Carlo (MC) simulations on Heisenberg 
spin systems have been carried out on a variety of spin 
lattices to examine the influence of magnetic dipole inter¬ 
actions on the properties with either dipolar interactions 
only or in combination with other spin interactions. For 
purely dipole interactions, Bouchaud and Zerah studied 
FM on the fee lattice and determined the Curie temper¬ 
ature Xc . 23 They studied the critical exponents at Tq 
and determined the anisotropy constants K\ and K 2 , 
where they found that the ordered moment direction in 
the collinear FM state at the lowest T was along [100], 
with a crossover from [111] at higher T. Tomita reported 
MC simulations on 2D triangular, square, honeycomb 
and kagome spin lattices with only dipolar interactions 
and studied the ground state magnetic structures and 
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critical phenomena. 24 One result was that the kagome 
lattice has a FM ground state with 1/3 of the spins dis¬ 
ordered at T = 0 [an amplitude-modulated (AM) state] 
with residual entropy (“missing entropy”) at T = 0 re¬ 
sulting from macroscopic degeneracy of the ground state. 
[An AM magnetic structure is one where the magnitude 
of the ordered moment is not the same for all (identical) 
spins in the spin lattice.] Very recent MC simulations on 
the kagome lattice by Holden et al. 25 and Maksymendo et 
al. 2e instead found a noncollinear coplanar equal-moment 
ground-state magnetic structure on the kagome latice. 
Thus when an AM magnetic structure is obtained the¬ 
oretically for a particular spin lattice, this may indicate 
that a lower-energy equal-moment magnetic structure ex¬ 
ists in which the moments have their maximum (satura¬ 
tion) value. 26 Other MC simulations examined the influ¬ 
ence of dipolar interactions on the properties in combi¬ 
nation with other interactions. 26-28 

In this paper previous work on the effects of mag¬ 
netic dipolar interactions on the magnetic and thermal 
properties of magnetic systems is significantly extended. 
Usually exchange and/or RKKY interactions are stronger 
than dipole interactions and determine the nature (FM or 
AFM) and k of the magnetic structure. However, when 
the exchange interactions are Heisenberglike (isotropic), 
some sort of anisotropy is needed to determine the di¬ 
rections of the ordered moments in the ordered state as 
discussed above, even if very weak compared to the ex¬ 
change interactions. The present work was initially mo¬ 
tivated by the lack of systematic studies of this topic for 
uniaxial tetragonal and hexagonal Bravais spin lattices 
versus the c/a ratio to compare with experimental re¬ 
sults such as for the Eu +2 spins S = 7/2 in EuCu 2 Sb 2 
on a bet sublattice that exhibit collinear AFM ordering 
below T N = 5.1 K. 29 ’ 30 

We study the influence of dipolar interactions on the 
magnetic ordering temperature T m , on the collinear or¬ 
dered moment directions and the temperature T depen¬ 
dence of the ordered moment and other properties at 
T < T m , and on the Weiss temperature in the Curie- 
Weiss law at T > T m in a systematic way for a variety 
of spin lattices including ID, 2D and 3D spin lattices us¬ 
ing our recent formulation of the Weiss MFT. 11 All spins 
in a given system are assumed to be identical and crys- 
tallograpIdeally equivalent. The 3D spin lattices studied 
here include sc, bcc, fee, simple tetragonal, bet and sim¬ 
ple hexagonal (triangular) Bravais lattices. Non-Bravais 
spin lattices are also studied which include the honey¬ 
comb (chickenwire) lattice, the kagome lattice and the 
Shastry-Sutherland lattice. For the uniaxial stacked lat¬ 
tices, the eigenvalues and eigenvectors of the magnetic 
dipole interaction (MDI) tensor are calculated for c/a ra¬ 
tios from 0.5 to 3. We utilize an appropriately modified 
theory to calculate the properties of noncollinear AFM 
structures and compare the results with calculations as¬ 
suming collinear AFM structures for the same k. Within 
MFT, the contributions of different sources of molecular 
fields to the Weiss temperatures and the magnetic order¬ 


ing temperatures are additive. Therefore, for example, 
when dipolar and exchange interactions are simultane¬ 
ously present, one can calculate the dipole contributions 
to good accuracy and then subtract them from the ob¬ 
served values to obtain the contributions from the ex¬ 
change interactions. Then with a model for the exchange 
interactions one can estimate their values. In addition to 
calculating the magnetic and thermal properties of pure 
dipolar magnets, the anisotropy in the susceptibility of 
Heisenberg AFMs in the PM regime with T > T m is also 
computed. We compare our predictions to the magnetic 
properties measured for illustrative real materials. In this 
paper we do not consider critical phenomena, domain 
formation and similar effects in FMs or other potential 
sources of magnetic anisotropy in a spin system such as 
single-ion effects. 

Our theoretical framework allows easy extensions to 
calculate the dipolar contributions to the magnetic prop¬ 
erties of spin lattices not discussed here such as collinear 
or noncollinear ordering on orthorhombic, monoclinic 
and triclinic Bravais or other non-Bravais spin lattices. 

In Sec. II we first write down the expressions relat¬ 
ing the macroscopic magnetic induction, applied mag¬ 
netic held and magnetization including shape (demag¬ 
netizing) effects. The expression for the local held seen 
by a spin is discussed in Sec. IIB. The part of that lo¬ 
cal held (the near held) due to discrete moments inside 
a macroscopic Lorentz sphere is discussed in Sec. IIC, 
together with the energy of a spin interacting with the 
near held. Applications of the general theory in Sec. IIC 
to magnetically-ordered states in collinear magnets, non- 
Bravais spin lattices and coplanar noncollinear helical or 
cycloidal AFMs are presented in Secs. IID, HE and IIF, 
respectively. The expression for the near held due to mo¬ 
ments within a Lorentz line (ID), circle (2D) or sphere 
(3D) is discussed in Sec. IIG. Some details about cal¬ 
culations of the MDI tensor are given in Sec. IIH. In 
Appendix A some information useful for implementing 
the theory in Sec. II is discussed. 

The calculations of the eigenvalues and eigenvectors 
of the MDI tensor for collinear magnetic structures with 
specihc magnetic propagation vectors for ID and 2D spin 
lattices are given in Secs. Ill, where the 2D spin lattices 
include the square, triangular, honeycomb and kagome 
lattices. Three-dimensional spin lattices are considered 
in Sec. IV, where results are given for the three cubic 
Bravais lattices, the two tetragonal Bravais lattices, the 
simple hexagonal lattice and the honeyomb lattice. For 
the 3D tetragonal and hexagonal lattices the eigenvalues 
and eigenvectors are obtained versus the c/a ratio from 
c/a = 0.5 to 3 in 0.1 increments. For all spin lattices, we 
carry out calculations of the MDI tensor of a central spin 
with its neighbors by direct summation with increasing 
radius away from the central moment until convergence 
is achieved within at least 0.001%. The convergence of 
the dipolar sums is discussed in the corresponding sec¬ 
tion, and representative convergence plots are given in 
Appendix B. 
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The predictions of the easy axis for collinear AFM or¬ 
dering are compared with experimental results for the 
simple-tetragonal Mn and Fe sublattices in BaMn 2 As 2 
and BaFe 2 As 2 , and for the bet spin lattices in GdCu 2 Si 2 , 
EuCu 2 Sb 2 and MnF 2 . For these cases we compare the 
results of the eigenvalues and eigenvectors versus the c/a 
ratio in graphical format with the experimental data, and 
the graphical results for other cases are placed in Ap¬ 
pendix C. The treatment of noncollinear AFMs is pre¬ 
sented in Sec. V, with application to the 120° ordering 
on the triangular lattice, to the 90° ordering on the dis¬ 
torted Shastry-Sutherland GdB 4 compound and to the 
undistorted 2D and 3D Shastry-Sutherland lattices. 

Section VI presents the calculation of the FM order¬ 
ing temperature Tq and AFM ordering temperature Tn 
arising from dipolar interactions within our recent for¬ 
mulation of MFT . 11 A quantitative discussion of the 
competition between FM and AFM ordering on cubic 
Bravais lattices versus the demagnetization factor of a 
sample in the absence of FM domain formation is given 
in Sec. VII. The properties of dipolar magnets in the 
magnetically-ordered state are derived in Sec. VIII. The 
ordered moment and heat capacity of dipolar magnets in 
zero magnetic field versus temperature are presented in 
Sec. VIII A, where the results are the same within MFT 
for both FMs and AFMs. The dipolar anisotropy param¬ 
eter K\ for uniaxial dipolar AFMs versus temperature 
is derived in Sec. VIIIB. Calculations of the perpen¬ 
dicular susceptibility below Tn and the associated criti¬ 
cal field for uniaxial AFMs are presented in Secs. VIIIC 
and VIIID, respectively. 

The Curie-Weiss law for dipolar magnets in the PM 
state is derived in Sec. IX, where the Weiss temperature 
is found to be anisotropic in general. In Sec. X we spe¬ 
cialize to spherical samples of collinear AFMs, where the 
anisotropic susceptibilities x f° r temperatures above Tn 
as well as both the parallel and perpendicular suscepti¬ 
bilities below Tn are presented and discussed. Examples 
of these anisotropic x(T) behaviors are given in Sec. X D 
for simple tetragonal lattices with c/a < 1 , c/a = 1 (sc) 
and c/a > 1 . 

The anisotropic %(T) of a Heisenberg-exchange AFM 
at T > Tn due to MDIs is derived in Sec. XI and applied 
to fit the experimental data for single-crystal MnF 2 . The 
paper concludes with a short summary in Sec. XII. 

Tables of values of the dipolar eigenvalues and eigen¬ 
vectors versus c/a plotted in the text and Appendix C 
are available in the Supplemental Material . 31 


II. THEORY 

The magnetization per unit volume of magnetic ma¬ 
terials can be significant compared to the applied field 
and results in a demagnetizing field and an internal field 
smaller than the applied field. In the following the theory 
for this important demagnetizing correction is discussed 
within the Gaussian cgs system of units 32,33 that is used 


throughout this paper. 

A. Macroscopic Fields 

We initially assume that a sample has the shape of 
an ellipsoid of revolution and that the applied field is 
along one of the three principal axes a. Then the volume 
magnetization (net magnetic moment per unit volume) 
M (units: G) is uniform in the sample and the magnetic 
induction B (units: G), the magnetic field H (units: Oe 
= G) and M are collinear with components M a , H a and 
B a for the external field H Q applied along the a axis. 
For each point in space one has 

B a = H a + 4nM a . (la) 

Thus internal to the sample one has 

-^int a. — a +47tM q . (lb) 

The demagnetizing field internal to the sample due to 
M a is 

H da = -4nN da M a , (lc) 

where here the demagnetizing factor N da is defined as 
in the Systeme International system of units for which 
0 < N da < 1 and -^dc* = 1- Thus the internal 

magnetic field H lnta and the magnetic induction B^ n ^ e 
due to sample shape effects and including the applied 
field H a are 

Hint a = H a ATrN da M a: (2a) 

S f„ h tT = H in t a +AnM a 

= H a + (1 - N da )4nM a . (2b) 

For a given M a , the internal field is Hi nta in Eq. (2a). 
Thus in descriptions of the magnetic behavior of a sample 
in terms of M a and H ai one can correct for the demag¬ 
netizing field by retaining the measured value of M a but 
replacing H a by H a — ATrN da M a , where N da is estimated 
from the sample shape and the field orientation with re¬ 
spect to the sample (see below). 

The magnetic susceptibility of a material is often 
defined as % = M(H)/H, which in general is field- 
dependent. In the present discussion, M is the volume 
magnetization, so x is the susceptibility per unit vol¬ 
ume and is dimensionless. The observed susceptibility 
is then x« bs = M a /H a and the intrisic susceptibility is 
Xa = M a /Hi n t a - Utilizing Eq. (2a) one obtains Xa from 
X° bs according to 

M a = M a = M a /H a 

Hmta H a — 4nN da M a 1 - 4-irN da M a /H a 

v obs 

— _An_ (%) 

1 - 47rV d aX° bs ‘ 

At each temperature one can correct the observed sus¬ 
ceptibility for the demagnetizing field using Eq. (3). 
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Alternatively, using Eq. (3) one can write the observed 
susceptibility in terms of the intrinsic one as 


X 


obs 

OL 


X* 

1 + A-KN Aa Xa' 


( 4 ) 


Thus when 4TrN Aa x a 1, one obtains the field- 
independent susceptibility and linear M a (H a ) behavior 


Y obs =_ 

4 7rJV da ’ 


M a 


4:irN da Ha ' 


( 5 ) 


The latter behavior holds until M a reaches it saturation 
(maximum) value M sat a ; at higher fields M a is of course 
equal to its constant saturation value. In practice, the 
limiting behaviors in Eqs. (5) are realized only when a 
material is approaching its FM transition temperature 
from above. 

An expression for the demagnetizing factors N da for 
the general ellipsoid of revolution was calculated long 
ago . 34 For sample shapes other than ellipsoids, M is 
not uniform within the sample except for limiting cases. 
What is then relevant in the present context is the de¬ 
magnetizing field averaged over the sample volume as 
expressed in the associated “magnetometric” demagne¬ 
tizing factor. Such sample shapes include the cylinder 
and the rectangular parallelepiped (rectangular prism) 
for which the magnetometric N Aa values have been calcu¬ 
lated for arbitrary sample dimensions in Refs. 35 and 36, 
respectively. 


B. Local Magnetic Induction from Magnetic 
Dipole Interactions 

Theoretical predictions of magnetic properties for local 
magnetic moments are often cast in terms of the local 
magnetic induction B{°“' seen by a local moment pi at 
position r^. This local magnetic induction along a given 
principal axis a is traditionally written for a 3D spin 
lattice in terms of the four contributions 37 

Ant + Ant a + Ant a + Ant ) 

where B a = H a is the applied magnetic induction aris¬ 
ing from currents outside the sample and Bf 4 ) 3 )) 6 is the 
contribution in Eq. (lb) due to the sample shape. The 
contribution 


Bf£T tz = yM a (7) 

is the Lorentz cavity field inside a spherical cavity of ra¬ 
dius R surrounding the point at its center at position r; 
at which is to be calculated. The fourth contribu¬ 

tion i is the sum of the dipolar fields at position 
arising from the other magnetic dipoles inside the Lorentz 
cavity at positions r ? . This is the only term that depends 
on the crystal structure of the material. The Lorentz 
cavity radius R is much larger than the distance between 
magnetic moments in a sample and is large enough so 


that the calculated f becomes independent of R to 
within some specified precision. Substituting Eqs. (2b) 
and (7) into (6) gives 

BlTai = H a + (i - N A ^j AirM a + B? n Zi. (8a) 

This is an important fundamental equation for calculat¬ 
ing the local field. 

Two special cases of Eq. (8a) are of use. In the first, one 
corrects the applied field for the demagnetizing field in 
the measurements as described above which is equivalent 
to removing N Aa from Eq. (8a), yielding 

A nr 

Blntai = H a + -M a + B%r ai . (8b) 

This equation is sometimes favored for comparison of 
theoretical predictions of the dipolar magnetic proper¬ 
ties with experimental data because it is independent of 
sample shape. Here M a is the total magnetic moment 
per unit volume. If all spins are identical and crystallo- 
graphically equivalent as assumed throughout this paper, 
one can write M a = fi a / I4 P i n where is the net average 
ordered and/or induced moment per spin in the a direc¬ 
tion and I4 P in is the volume per spin, so an equivalent 
form of Eq. (8b) is 

A -TT 

= H a + —— n a + B%? ai . (8c) 

•3 ^spin 

Note that M a = [i a = 0 for an AFM in H = 0. 

Alternatively, one can shape a sample into a sphere, 
giving N Aa = 1/3 for all three principal directions a, and 
then Eq. (8a) becomes 

plocal _ TT < pnear /'CU'l 

-“int ai ~ + ^int ail 

which eliminates the effect of the Lorentz field but only 
applies to a spherical sample. This formulation is de¬ 
sirable if one wishes to ameliorate the tendency of the 
Lorentz field to enhance dipolar FM ordering with re¬ 
spect to AFM ordering, as illustrated in Fig. 13 below 
where FM is favored for small values of N Aa for bcc and 
fee Bravais lattices. 


C. Magnetic Induction Due to Collinear Alignment 
of Magnetic Dipoles Inside Lorentz Cavity 

The magnetic induction B (J seen by a central moment 
pi at a position due to a point magnetic dipole moment 
pj at position r ? is 

B u = -g-[3 {Pj ■ v,, 'lb- - rfipj], (9a) 

r ji 


where 


r,; = r, - r, 


(9b) 
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The energy of interaction Ei of pi at position r,, due 
to the magnetic induction By is 

Ei = — —/Xj ' = 2 r 5 ' r J®) — ' AbL 

( 10 ) 


where the factor of 1/2 in the first equality recognizes 
that the interaction energy of the jli with By from /X, is 
equally shared between pi and pj. Expanding the first 
term on the right side of Eq. (10) in Cartesian coordi¬ 
nates, one can write the term in matrix form as 


J 


Tt 


ifl'i * C i ) (Ab ' U i ) — (/O'.T Aby Ab 2 ) I Tiicr.T 

1 jix 1 jiz 


jix 

jix r jiy r jiy 

T* irxX ii.7. T i ii/V 


jix'jiy 1 jix 1 jiz 


T iiii'T ii 


jiy 1 jiz 

jiy jiz r ji z 


fJ’jx 

A 4 jy | : A 4 / r ji r jiPji 

Hjz 


( 11 ) 


where pj is the transpose of the column vector pi, pj is 
a column vector and is a 3 x 3 dyadic. Similarly, 

the scalar product in the second term on the right side 
of Eq. (10) can be written in matrix form as 

A L i ' Ab = Ah lAb> (12) 


both FM (k = 0) and AFM structures below the magnetic 
ordering temperature T m and the FM-aligned magnetic 
structure induced above T m by an external magnetic field 
applied along one of the three principal axes of the MDI 
in Eq. (16c) below. From Eq. (14) one obtains the “ex¬ 
tinction condition” 


where 1 is the 3x3 identity matrix. Using Eqs. (11) 
and (12), Eq. (10) can be summed over all neighbors pj 
within a length of chain (ID), a circle of specified radius 
(2D) or Lorentz sphere (3D), all centered on pi, and then 
can be succinctly written in matrix form as 

Ei = -^pjGipj. (13a) 

where the 3x3 symmetric tensor Gi is 

Gi = J2 i( 3r A r A - rfi 1 )- ( 13b ) 

In order to solve Eq. (13a) for the eigenenergies Ei and 
eigenvectors pi of the tensor G,, one must first express 
each pj in terms of pi- In the following three sections we 
discuss our methods for doing so for collinear magnetic 
structures on Bravais and non-Bravais spin lattices and 
coplanar noncollinear AFM structures, respectively. 


D. Collinear Magnetic Structures 

In this section we consider collinear magnetic struc¬ 
tures with magnetic wavevector k where 

pj = cos(k • r ji)pi. (14) 

Since the cosine function is a scalar with a value between 
±1, Eq. (14) expresses that pj can be either parallel or 
antiparallel to pi. For cos(k • TjP) = ±1 for each pj the 
magnetic structure is an “equal-moment” (EM) structure 
where the ordered moments all have the same magnitude 
/x (which depends on T). For cos(k ■ r 7 , : ) ^ ±1 for some 
pj, [i depends on j and the structure is a collinear AM 
AFM structure. Collinear magnetic structures include 


Pj — 0 if k • r ji = odd multiple of — rad, 


(15) 


as in Eq. (A 6 ) for AM AFM structures associated with 
specific k values and spin lattices. A general k corre¬ 
sponds to either an EM or AM collinear AFM structure. 

All simple Bravais lattices have EM magnetic struc¬ 
tures. Amplitude-modulated AFM structures occur 
when the simple Bravais lattices have more than one 
spin in the unit cell such as for bcc, bet and fee spin 
lattices. With the cos(k • r ji) term as given in Eq. (14), 
EM structures occur for bcc and bet lattices with k = 
(i,i,0), (i,0, i) and (001) r.l.u., and AM structures 
for k = (A, 0,0) and (|, |, i) r.l.u. For the fee lattice, 
EM structures occur for k = (i, ;, -■ ) and (0,0,1) r.l.u, 
whereas AM structures occur for k = (A, 0,0), (^, 0), 

(i,0, l) and (|, |, i). With the exception of the last one, 
all AM structures considered can be converted into EM 
structures by inserting an additive phase in the cosine 
term: cos(k ■ r ji) —► cos(k ■ r ji + </), where (j) = 7t/4 rad. 
In that case, all eigenvalues are reduced in magnitude by 
the factor cos(7t/4) = 1 /\/2, which corresponds to a re¬ 
duction in the ordered moment by a factor of 1/2 1 / 4 . All 
eigenvalues plotted or listed in this paper were obtained 
for cj) = 0. 

In pure magnetic dipole AFMs, the above discussion 
shows that the AFM ground state can be an AM state, 
depending on the AFM wavevector. However, even in 
systems in which the magnetic dipole interaction is not 
expected to play an important role, this interaction can 
still cause a small modulation of the ordered moment 
versus position in the magnetic unit cell. Furthermore, 
large-amplitude AM AFM structures are observed in 
geometrically-frustrated systems such as in Gd 2 Ti 2 07. 38 
Because AM structures contain at least some fraction of 
spins with ordered moments less than the saturation mo- 
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ment and hence show strong quantum fluctuations in the 
ground state, the entropy increase on heating from low 
temperatures would be less than the value R 111(25 + 1) 
per mole of spins. This can be checked by calorimetry. 

The discussion throughout this paper applies to iden¬ 
tical crystallographically-equivalent spins with identical 
saturation moments /r sat and with thermal-average (or¬ 
dered) magnetic moments /q = [ifa, where p can be dif¬ 
ferent for different spins in AM structures. We express 
rji in units of the lattice parameter a of the respective 
crystal structure. The crystallographic unit cell often 
contains more than one spin per unit cell in the exam¬ 
ples described. Then using Eq. (14), Eqs. (13) become 39 


that the MDI at least contributes to ordering along the 
observed easy axis; however, this does not rule out other 
sources of anisotropy that may also contribute. 

A general feature of the eigenvalues Ak Q of the MDI 
tensor G, for a given k and spin lattice is that their sum 
over the three eigenvectors a is identically zero when no a 
priori constraint is placed on the ordering axis of /q. This 
sum rule is violated when such a constraint is imposed 
such as for coplanar noncollinear helical or cycloidal AFM 
order as discussed in Secs. IIF and V A. In those cases, 
one of the Aka is the eigenvalue for FM ordering (k = 0) 
along the helix or cycloid axis. The other two eigenvalues 
and corresponding eigenvectors are the ones associated 
with the actual AFM components of the helix or cycloid. 


Ei = -e (16a) 

where 


E. Non-Bravais Spin Lattices 


e = 



(16b) 


has dimensions of energy and the dimensionless symmet¬ 
ric MDI tensor is 


Gi ^ ( W «) 5 



cos(k • r^). 


(16c) 


Labeling the eigenvalues of Gj(k) as Aka, Eq. (16a) gives 
the eigenenergies as 


Ei a — e Aka, (16d) 

where the subscript a refers to a Cartesian principal or¬ 
dering axis eigenvector of the collinear magnetic struc¬ 
ture, where the three principal axes are orthogonal to 
each other. Thus the ground state energy and order¬ 
ing axis for a given k due to the MDI corresponds to 
the largest of the three Aka eigenvalues. The MDI en¬ 
ergy scale is set by the value of e in Eq. (16b) which is 
system-dependent. The value of e/fce is typically of order 
0.01-0.1 K. 

The magnetic propagation vector k must be specified 
in terms of the reciprocal lattice translation vectors in 
Cartesian coordinates in advance of computing Gj a (k). 
One can calculate the Aka eigenvalues and correspond¬ 
ing eigenvectors (ordered moment axes /q) for various k 
vectors, including k = 0 for FM-aligned moments which 
may occur due to FM ordering in applied field H = 0 or 
to H > 0 in the PM state. Usually the magnetic k vec¬ 
tor observed by, e.g., neutron diffraction measurements, 
is determined by exchange or RKKY interactions rather 
than dipole interactions. In that case one can still test 
whether the easy axis predicted by the MDI is consistent 
with the observed one. A negative answer would indi¬ 
cate that the MDI does not contribute to determining 
the easy axis, and hence some stronger source of magne¬ 
tocrystalline anisotropy must be present that overcomes 
the preference of MDIs. A positive answer would mean 


A crystal structure consists of a Bravais lattice plus 
a basis of atoms attached to each Bravais lattice point. 
Non-Bravais spin lattices are Bravais lattices with more 
than one spin in the basis. These include, e.g., the fee 
diamond lattice and the 2D hexagonal honeycomb (or 
chickenwire) lattice, each with two spins in the basis, 
and the kagome lattice with three spins in the basis. In 
such cases one must modify Eq. (16c) to include a sum 
over the atoms in the basis, in addition to the sum over 
Bravais lattice points already included in Eq. (16c) via, 
e.g., Eqs. (Al). AFM structures in such non-Bravais spin 
lattices include those with AFM propagation vector k = 
0 for Neel-type ordering on the 2D honeycomb lattice, 
which is the same propagation vector as for FM order¬ 
ing. In such AFM structures where the magnetic and 
crystallographic unit cells are the same, in order to cal¬ 
culate Gi one must specify the orientations of the ordered 
moments within a unit cell with respect to the orienta¬ 
tion of a central moment /I). Thus Eq. (16c) is modified 
to read 




(rjki/a) 


, YjkiYjki 
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(17a) 


where the sum over j again refers to the sum over the 
Bravais lattice positions, the sum over k sums over all 
atoms in the basis, the position of the central moment 
jli is not necessarily at the origin of of a central unit cell, 
and the vector from jli to a moment jlk is 


r jki = +rfc - rj, (17b) 


where is the position of moment jlk in the basis with 
respect to the position of the associated Bravais lattice 
point Tj. The term with rjki = 0 is omitted from the sum 
because that term corresponds the difference in position 
of moment jli with itself. The Cartesian rotation matrix 
Rki in Eq. (17a) expresses the moment direction of jlk in 
the basis with respect to that of the central moment jli 
via 


k'k — 


(17c) 






similar to Eq. (14) for collinear ordering associated with 
a magnetic propagation vector k. Prior to calculating 
G,, the 3x3 Rfci rotation matrix must be specified for 
each spin in the basis via a model for the AFM structure. 
For example, for the Neel AFM structure in Fig. 1 below, 
if fli were at a red position r i/a = |a+ |b, then R|for 
a spin at another red position would be R/^ = 1 and 
that for a black position would be —1, where again 1 
is the 3x3 identity matrix. This procedure is easily 
generalized to more than two spins per Bravais lattice 
point, as illustrated in Sec. VB below for calculating G; 
for the known coplanar noncollinear AFM structure of 
tetragonal GdB 4 in Fig. 12 containing four moments in 
the basis, each pointing in different directions, and for 
the related Shastry-Sutherland spin lattice. 


F. Coplanar Noncollinear Helical or Cycloidal 
Antiferromagnets 

Here we extend the above discussion to coplanar non¬ 
collinear helical or cycloidal AFM ordering on tetragonal 
or hexagonal Bravais lattices. For both types of AFM 
order, the ordered moments are defined to he in the crys¬ 
tallographic ab plane. For helical AFM ordering, the or¬ 
dered moments are FM-aligned in the ab plane and the 
helix wavevector k axis is the c axis. For cycloidal AFM 
ordering, k lies in the ab plane and the moments in planes 
perpendicular to both k and the ab plane are aligned fer- 
romagnetic.ally. The Cartesian ar-axis is parallel to a, 
the y-axis is perpendicular to a in the ab plane and the 
z axis is perpendicular to the ab plane along the c axis. 
Pictures of the helical and cycloidal structures are given 
in Refs. 40 and 41, respectively. In either structure, the 
azimuthal angle 4 >ji = 4>j — (f>i with respect to the posi¬ 
tive a axis between moments /2 ? and pi in the ab plane 
is given by 


Pji — k • Tji. 


(18a) 


The relationship between the central moment direction 
jli at position r i and that of another moment at position 


i'j in either AFM structure is 


0 0 


N = 


0 sin (f>ji 0 | jli 

0 1 


COS (f>ji 

0 

which can be written 

pj = (xx cos (j)ji + yy sin fa + zz )/q, 


(18b) 


(18c) 


where the Cartesian coordinate system is used through¬ 
out. Then Gj in Eq. (16c) becomes 


G* = E , 

v 0 


fai/a ) 5 


r • r • - V ■ ■ 

o L 3i L ji _ 3* 

a 2 a 2 


(19) 


x (xx cos fa + yy sin^ + zz). 


As with collinear AFM ordering, one must specify k 
in terms of the reciprocal lattice translation vectors in 
Cartesian coordinates in advance of computing Gj(k). 
Note that when Gi(k) is diagonalized, one eigenvalue and 
corresponding eigenvector are for FM ordering along the 
z axis and are not relevant to those for the helix, whereas 
the other two sets of eigenvalues and eigenvectors are for 
the helix. As a result, the sum of the three eigenvalues do 
not add to zero as they do for all other AFM structures 
discussed above. 


G. Near Field 


The value of B^C in Eq- (6) that is seen by a given 
moment pi at position r, in a given magnetic structure 
with a given ordered moment configuration, due to the 
sum of the magnetic fields from the magnetic moments 
around it within the Lorentz cavity of radius R , is simply 
given as 


B 


near 
int a i 


2Ei fl\ ko; 

fi a a 6 


( 20 ) 


where the factor of 2 arises because the energy per pair 
is split evenly between each pair of moments, whereas 
the magnetic field arises only from the neighbor of each 
pair, the second equality was obtained using Eqs. (16b) 
and (16d) and B-J^C can be either positive or negative, 
depending on the sign of Aka- If the MDI is the only 
source of anisotropy present, this field must be positive 
because then the ordered moment is parallel to the local 
magnetic induction, which minimizes the free energy of 
the moment. The quantity B?® a G is needed to calcu¬ 
late the total local magnetic induction at the site of a 
local moment according to Eq. (8a). If £) is expressed 
in cgs units of erg and those of fj, in cgs units of erg/G 
(= Gem 3 ), then B“ t a C has the correct cgs units of G. 

Using Eq. (20), the total local field in Eq. (8b) seen by 
central moment pi becomes 


zolocal _ TT 

n intai ~ 11 a 


A-Kflg/fl 

3Uspin/a 3 


+ Aka ) —o) (21) 


where the first term in parentheses is the Lorentz field, 
where we distinguish the moment component fi a in the 
a direction per spin averaged over the sample and the 
magnitude fj, of the average moment per spin. A nonzero 
value of fi a only occurs in a FM or in an AFM in the 
presence of an external magnetic field. The second term 
in parentheses arises from the near held. 


H. Calculation and Diagonalization of the 
Magnetic Dipole Interaction Tensor 

We chose to carry out the sums in the expressions 
for the dipole interaction tensor G^ in Eqs. (16c), (17a) 
and (19) directly instead of by using the Ewald-Kornfeld 
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technique , 42 because we wanted to study the convergence 
properties of the eigenvalues Aka versus the radius R of 
the circle or Lorentz sphere for 2D and 3D lattices, re¬ 
spectively. The calculations and diagonalizations of G, 
were carried out using standard Macintosh laptop and 
desktop computers and Mathematica software. For the 
ID chain with FM and Neel AFM states, the eigenvalues 
and eigenvectors of Gj are trivially determined exactly 
for the infinite chain as shown in Sec. Ill A. For 2D lat¬ 
tices the sums were carried out within circles of radius 
up to R/a = 1000 containing up to 1 x 10 7 spins (for 
the kagome lattice containing three spins per unit cell). 
For the 3D lattices the sums were carried out within a 
Lorentz sphere, usually up to a radius R/a = 50 contain¬ 
ing up to 6 x 10 6 spins. Calculations were also done for 
two AFM structures out to a sphere radius of R/a = 100 
containing 1.7 x 10' spins to check convergence. 

For the FM spin structures in 2D, the values of Gj 
versus l/(R/a) were extrapolated to 1 /(R/a) = 0. As 
shown in Appendix B, the calculations of Gj for AFM 
structures generally converge more rapidly with increas¬ 
ing R/a than for FM structures. These procedures deter¬ 
mined Aka to accuracies of < ±10 ~ 6 for 2D lattices and 
< ±0.001 for 3D lattices, more than sufficient for our 
purposes. The eigenvectors jli usually converged very 
quickly with increasing R/a. For the various 3D tetrag¬ 
onal and hexagonal lattices, Gj was calculated for c/a 
ratios from 0.5 to 3 in 0.1 increments. 

Figures 17 and 18 in Appendix B show the convergence 
of Aka with increasing R/a for FM and Neel AFM mo¬ 
ment alignments along the c axis in the 2D simple square 
lattice, respectively. Figures 19(a) and 19(b) show plots 
for a simple tetragonal lattice with FM alignment of the 
moments along the c axis for c/a = 1.5 and 3, respec¬ 
tively. Figures 20(a) and 20(b) show analogous plots for 
the simple tetragonal lattice with Neel-type AFM order¬ 
ing where k = (§, b , I) and alignment of the moments 
along the c axis for c/a = 1.5 and 3, respectively. 

For the 3D FM and AFM structures, the values of Aka 
were typically obtained for R/a = 1 to 50 in increments 
of 1 and the last 10 or 20 values were averaged to obtain 
the data in the figures in the text and Appendix C and 
in the tables in the Supplemental Material . 31 


III. EIGENVALUES AND EIGENVECTORS 
FOR MAGNETIC ORDERING ON ONE- AND 
TWO-DIMENSIONAL SPIN LATTICES 


A. Spin Chain 


We assume that the spin chain lattice is oriented along 
an axis desigated as the a axis (x axis) with spacing a 
between adjacent spins, so 

TJL = n a . ( 22 ) 


Ferromagnetic alignment corresponds to k = 0. This 
alignment can occur either in the FM-ordered state or 
in the PM state in the presence of an applied mag¬ 
netic field. The central spin is positioned at n a = 0, 
so n a of the neighbors runs from —oo to oo, excluding 
n a = 0. Numerical diagonalization of Gj in Eq. (16c) 
with k = 0 and |n™ ax | = 1000 (2000 neighbors of the 
central moment) shows that the principal axes of the 
interaction tensor are parallel and perpendicular to the 
a axis. The lowest energy configuration with a calcu¬ 
lated A(o,o,o)[iool = 4.80823 is with the ordered moments 
aligned along the a (chain) axis. This makes sense be¬ 
cause the lowest energy configuration of a moment is 
when each moment points along the local field seen by 
the moment, which is along the axis of the chain. The 
eigenvalues for the two higher-energy orthogonal direc¬ 
tions are A( 0 ,o,o)[oio] = -\o,o,o)[ooi] = — \o,o,o)[ioo]/2- 
For the present case of the ID spin chain one can also 
evaluate A( 0 .o,o)[i.o,o] exactly. Equation (16c) yields the 
eigenvalue 


A(0,0,0)[1,0,0] = 4 Z3- ( 23 ) 

n a — l ' a 

The sum is = C(3 ), 43 yielding 

A ( o,o,o)[i,o,o] =4C(3) w 4.808 228 (24) 

as shown in Table I, where £(z) is the Riemann zeta 
function with £(3) ss 1.20206. The above numerical value 
of 4.80823 obtained for A( 0 ,o.o)[i,o.o] agrees with this exact 
value to six-place accuracy. This shows that the value of 
|n“ ax | = 1000 and a spin chain containing 2000 neighbors 
of the central moment used in the numerical calculation 
is sufficient to obtain this accuracy. 

It is of interest to examine the approach to the infinite- 
chain limit of A(o,o,o)[i.o,o] on increasing \n a \- For large 
\n a \ one can replace the sum in Eq. (23) in the region 
where n a is large by an integral f n~ 3 dn a oc —1/n 2 . Thus 
we expect for n a 1 that 


A(o,o,o)[i,o,o] — 4C(3) 


A 


4C(3) 


A ' 
403KL 


(25) 


where A is a positive constant. An exact series ex¬ 
pansion of the sum in Eq. (23) about n a = oo indeed 
gives A( 0 ,o,o)[i,o,o] = 4£(3) - 2 /n\ + 0 (n~ 3 ), yielding 
A = 2. Equation (25) then predicts six-place accuracy 
for A(o t o,o)[i,o,o] f° r \ n a\ = 1000, consistent with the above 
comparison. 

Here we also examine the Neel-type AFM wavevector 
k = ( 1 / 2 , 0 , 0 ) r.l.u., where 1 r.l.u. = 2 n/a is the re¬ 
ciprocal lattice unit for this spin lattice. A numerical 
calculation using Eq. (16c) shows that the eigenvalues of 
Eq. (16d) converge to six significant figures even with a 
small |n a | max = 70. These calculations also show that 
the most stable ordered moment direction is perpendic¬ 
ular to the chain with 


A(i/2,o,o)[o,i,o] — A(i/2,o,o)[o,o,i] — 1.80309 (26) 
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and the unstable x-axis direction has 


•^( 1 / 2 , 0 ,o)[i,o,o] — ~ 2 A(i/ 2 ,o,o)[o,i,o] — —3.60617. (27) 

An exact calculation for the a-axis eigenvalue is ob¬ 
tained using the AFM version of Eq. (23), yielding 


^(1/2,0, o)[i, o,o] 



-3C(3) « -3.60617, 
(28a) 


as listed in Table I. This value is identical to within 
six places with the numerical result for A(i/2,o,o)[i,o,o] ob¬ 
tained above using only a 141-spin chain (including the 
central spin). Thus the dipole fields seen by a central 
moment converge much faster with increasing n™ ax for 
the AFM structure than for FM one for the spin chain. 
The two lower-energy eigenvalues are 


^( 1 / 2 , 0 , 0 ) 10 , 1 , 0 ] 


^(l/ 2 , 0 , 0 )[ 0 , 0 ,l] — ~ 2 ^( 172 ,o, 0 )! 1 ,°,o] 
|C(3) « 1-80309. (28b) 


Comparing Eqs. (24) and (28b) shows that the eigen¬ 
value for the FM-aligned state with the ordered/induced 
moments aligned along the a axis is larger than the max¬ 
imum AFM one, and hence the energy per spin is smaller 
according to Eq. (16d) for the FM state than for the AFM 
state. The FM state is thus expected to be the magnetic 
ground state of the linear spin chain for purely dipolar 
interactions provided the ordering is not destroyed by 
quantum fluctuations. 


B. Two-Dimensional Spin Lattices 

For the simple square lattice, one has a = b and the 
spin positions given by the first two terms in Eq. (Ala). 
The normalized wavevectors are the first two terms in 
Eq. (A4). The largest (positive) eigenvector Aka of 
Gja(k) in Eq. (16c) with k = 0 for FM alignment 
occurs for the a or b easy axes, with A( 0 ,o,o)[ 1 , 0 , 0 ] = 
A(o,o,o)[o,i,o] an d A(o,o,o)[o,o,i] = —2A(o,o,o)[ 1 , 0 , 0 ]- Shown 
in Fig. 17(a) in Appendix B is the dependence of 
A( 0 ,o,o)[o,o.i] on the inverse radius R~ l for the c-axis 
eigenvalue A(o,o.o)[o,o,i]- According to the discussion in 
Sec. Ill A, in 2D one should have A( 0 ,o,o)[o,o,i](-R/ a ^ 
1) = A( 0 ,o.o)[ 0 , 0 , 1 ] ( a /= 0 ) + A/R, in agreement with 
the calculations in Fig. 17(a). Fitting the data for 
0.001 < a/R < 0.002 gives 

A(o,o,o)[o,o,i](a/i? = 0) = -9.0336220(1), (29) 

A = 6.28356(9). 

The deviations of the data from the fit are shown in 
Fig. 17(b), where it is seen that the deviations are of 
the order of 1 part in 10' for 0.001 < a/R < 0.0026. The 
graininess of the lattice becomes more apparent at larger 
values of a/R. 



FIG. 1: (Color online) Honeycomb lattice. Each honeycomb 
cell (not a unit cell) is bounded by solid blue lines. The hexag¬ 
onal unit cells with translation vectors a and b are outlined 
by dashed black lines. The 2D space group is p6m (No. 17) 
with two spins in Wyckoff positions 2b (§,§), (|,|). Bipar¬ 
tite Neel ordering of the two spins per unit cell is shown. The 
solid red circles represent half the magnetic moments pointing 
in one direction and the open black circles correspond to half 
the moments pointing in the opposite direction. 


The eigenvalues and eigenvectors for square-lattice 
AFM propagation vectors k = (d, 0, 0) (stripe AFM) and 
(|, 0) (Neel-type AFM) were also computed as shown 

in Table I. One sees that of these and the FM propaga¬ 
tion vector, the stripe AFM propagation vector with the 
ordered moments aligned along the b axis (perpendicular 
to k as shown in the last column of Table I) has the low¬ 
est energy. Our eigenvalues A ( 0 ,o,o)[ioo] and A ( 1 0 ,o)[ioo] 
are in agreement with, but are more precise than, those 
previously reported in Ref. 44, and our A(o,o,o)[ooi] an d 
A(i/ 2 , 1 / 2 , 0 )[ 001 ] values are in precise agreement with the 
results in Ref. 45. 

For the 2D simple-hexagonal (triangular) lattice the 
eigenvalues and eigenvectors were calculated for the FM 
state and four AFM states. From Table I, the lowest- 
energy states are the FM state and the AFM state with k 
= ( 1 , 0 , 0 ) (stripe-type), with the moments aligned within 
the ab plane in both cases. Data for the AM AFM state 
with k = (i, i, 0 ) are included in Table I because the 
classical ground state of a triangular lattice of spins with 
Heisenberg interactions is the well-known coplanar non- 
collinear 120 ° state that can be described by k = (|, 0 ) 

which we consider further in Sec. V A. Our eigenvalues 
A(o,o,o)[ioo] an d A^ 1 i 0 )[_ 1 / 2 l - A / 3 / 2 , 0 ] are in agreement to 
seven significant figures with those previously reported 
in Ref. 46. 

The 2D hexagonal honeycomb lattice is a non-Bravais 
spin lattice containing two spins per unit cell as shown 
in Fig. 1. The eigenvalues and eigenvectors of G; for 
collinear magnetic ordering were calculated for this lat¬ 
tice according to the method of Sec. IIE and the results 
are listed in Table I for the FM and Neel-type (Fig. 1) 
AFM states, where the magnetic propagation vector for 
both states is k = ( 0 , 0 , 0 ). 

Rozenbaum found that the ground state of the 2D 
honeycomb lattice is noncollinear with all spins aligned 









11 


TABLE I: One- and two-dimensional spin lattices. Eigenvalues Ak« and eigenvectors /t = [fi x , fj, y , y, z \ in Cartesian coordinates 
of the MDI tensor G;(k) in Eq. (16c) for various values of the magnetic wavevector k in reciprocal lattice units with collinear 
magnetic moment alignments. The most positive Aka value(s) corresponds to the lowest energy value according to Eq. (16d). 
The Cartesian x, y and z axes are along the a, b and c axes of orthogonal-axis lattices, respectively. For the hexagonal lattice, 
the x axis is parallel to the hexagonal a axis and the y axis is perpendicular the the a axis, rather than along the b axis. The 
linear chain is aligned along the a axis and the square and hexagonal lattices are aligned in the ab plane. 


k 

Ak/i 

A 

fi ■ k 

ID linear chain 

(0, 0, 0) (FM) 

4C(3) « 4.808 228 
-2C(3) « —2.404 114 

[100] 

[010], [001] 


(i, 0, 0) (Neel AFM) 

§C(3) w 1.803 085 

[010], [001] 

0 

-3C(3) « -3.606 171 

[100] 

1 

2D square lattice 

(A, 0, 0) (stripe AFM) 

5.098 873 

[010] 

0 

0.935 462 

[001] 

0 


-6.034 335 

[100] 

1 

(0, 0, 0) (FM) 

4.516 811 
-9.033 622 

[100], [010] 

[001] 


(i, i, 0) (Neel AFM) 

2.645 887 

[001] 

0 


-1.322 943 

[100], [010] 

l/%/2« 0.7071 

2D simple hexagonal (triangular) lattice 

(0, 0, 0) (FM) 

5.517 088 
-11.034176 

[100], [010] 

[001] 


(1, 0 , 0) 

5.517 088 

[100], [010] 

1 


-11.034176 

[001] 

0 

(|, 0, 0) 

4.094 909 

[i-#,0] 

1/2 


1.839 029 

[001] 

0 


-5.933 939 

0] 

—\/3/2 « -0.8660 

(i i 0) 

4.094 909 

[-i-^.o] 

— (b3+ l)/2 3/2 « -0.9659 


1.839 029 

[001] 

0 


-5.933 939 

[-0] 

l/b2« 0.7071 

(i i o) 

2.331 796 

[001] 

0 


-1.165 898 

[100], [010] 

l/b2« 0.7071 

2D hexagonal honeycomb lattice 

(0, 0, 0) (FM) 

17.092 359 
-34.184 718 

[100], [010] 

[001] 


(§, 0 , 0) 

12.827 051 

[-5.^.0] 

-1/2 


-0.090 183 

[001] 

0 


-12.736 868 

o] 

b3/2 « 0.8660 

(0, 0, 0) (Neel-type) 

12.116 366 
-6.058 183 

[001] 

[100], [010] 


2D hexagonal kagome lattice 

(0, 0, 0) (FM) 

51.321197 

11.205 800 
-62.526 996 

[010] 

[100] 

[001] 


(0, 1, 0) (ferrimagnet) 

40.458 644 

[001] 

0 


-0.171 624 

[100] 

0 


-40.287 021 

[010] 

1 

(§. I, o) 

13.213 509 

[001] 

0 


9.212 253 

[100] 

l/b2« 0.7071 


-22.425 762 

[010] 

1/n/2« 0.7071 

(0, i, 0) 

4.094 910 

[100] 

0 


1.839 029 

[001] 

0 


-5.933 939 

[010] 

1 
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AAAA 




FIG. 2: (Color online) Two-dimensional hexagonal kagome 
lattices. The hexagonal unit cell is shown at the lower left 
of each panel oulined in heavy black lines and contains three 
spins. The unit cell edges a and b are twice the length of the 
triangular-lattice unit cell edge. Three magnetic structures 
are shown, (a) The red, blue and open circles represent mo¬ 
ments that are mutually at an angle of 120° to each other 
within the ab plane, so a given moment has no nearest neigh¬ 
bors with the same orientation. The cycloid spin configura¬ 
tion shown has a wavevector k = (§, §,0) r.l.u. (b) Collinear 
magnetic structure for k = (0,1, 0) r.l.u. The red circles rep¬ 
resent moments in one direction and the blue circles repre¬ 
sent moments in the opposite direction. Because there are 
twice as many red as blue circles, this magnetic structure is 
a ferrimagnet (a net FM). (c) Collinear AFM structure for 
k = (0,1,0) r.l .u. The red and blue circles have the same 
meanings as in (b). There are equal numbers of red and blue 
circles, but the open circles represent spins with zero ordered 
moment, so the magnetic structure is an AM AFM. 


in the ab plane , 46 in contrast to the collinear FM and 
AFM structures assumed in the above calculations. He 
gave the ground state energy per spin as E/N = 
— 2 ^ 3 “ (4-453 809), where a nn = a/V 3 is the nearest- 
neighbor spin-spin distance. Converting to our notation 
for N = 2 spins per unit cell according to Eqs. (16) gives 
the eigenvalue A = 3 3 / 2 2(4.453809) = 46.2853. This 
eigenvalue is more than a factor of two larger (more sta¬ 
ble) than the most stable collinear magnetic structure in 
Table I for the 2D honeycomb lattice. 

The 2D hexagonal kagome lattice is very popular for 
studying the effects of geometric frustration on the prop¬ 
erties for AFM interactions. This lattice is shown in 
Fig. 2. The lattice is generated from a doubled trian¬ 
gular lattice by removing the spin at the origin of the 
unit cell, which shows that a kagome spin lattice is a 
^-depleted triangular lattice containing three spins per 
hexagonal unit cell. The 2D hexagonal space group of 
the kagome lattice is p6m (No. 17), with three spins in 
Wyckoff positions 3c (^,0), (0, |) , ( 5 , 5 )- For the 
kagome lattice the cycloid wavevector in Fig. 2(a) is k = 
(|, |, 0 ) r.l.u. instead of k = (A, |, 0 ) r.l.u. for the trian¬ 
gular lattice, due to the factor of two increase in the el¬ 
and 5-axis lattice parameters compared to the triangular 
lattice. The magnetic structure shown in the figure is the 
well-known classical 120 ° structure for nearest-neighbor 
AFM Heisenberg interactions. However, the ground state 
for collinear magnetic ordering arising from only dipole 
interactions is seen from Table I to be a FM structure 
with the moments pointing perpendicular to the plane of 
the lattice. Also shown in the table are results for two 
AFM wavevectors directed along the b* (y) direction dis¬ 
cussed next. 

A net FM (ferrimagnetic) collinear structure is shown 
in Fig. 2(b) with magnetic wavevector k = (0,1,0) r.l.u. 
There are twice as many moments pointing one way com¬ 
pared to the other way, so the net ordered FM moment is 
H S at/3, where // sa t is the saturation moment of each spin. 
An AM collinear AFM structure is shown in Fig. 2(c) 
with AFM propagation vector k = (0, A,0) r.l.u. The 
red and blue circles have the same meaning as in (b), but 
the black open circles represent spins with no ordered mo¬ 
ment. Therefore the average AFM ordered moment per 
spin is 2/i sat /3. 

According to Table I, the lowest-energy (largest eigen¬ 
value) collinear magnetic structure for the 2D kagome 
lattice is the FM structure with moments directed along 
the y direction (vertically upwards in Fig. 2) within the 
ab plane. The collinear structures shown in Figs. 2(a)- 
(c) are significantly less stable. Classical MC simula¬ 
tions determined that the ground state magnetic struc¬ 
ture is an EM noncollinear ferrimagnetic structure with 
all ordered moments lying in the ab plane . 25 ’ 26 The 
ground state energy per spin is quoted as E/spin = 
—2.38895 /i 2 /a 3 n , where a nn = a/2 and a nn is the 

nearest-neighbor spin-spin distance . 25 In terms of our no- 
2 

tation, E/spin = — //rr A which takes into account the 
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TABLE II: Simple cubic spin lattice. Eigenvalues Aka and 
eigenvectors fl = [fi x , y y , fi z ] in Cartesian coordinates of the 
MDI tensor Gj(k) in Eq. (16c) are listed for various values 
of the magnetic wavevector k in reciprocal lattice units with 
collinear magnetic moment alignments. The most positive 
Aka value(s) corresponds to the lowest energy value accord¬ 
ing to Eq. (16d). Also shown are the differences between the 
eigenvalues for the different eigenvectors for a given k and 
spin lattice, which are proportional to the respective mag¬ 
netic anisotropy energies and fields. The Cartesian x, y and 
z axes are along the a, b and c axes of the cubic unit cell, re¬ 
spectively. The labels A-, B-, C- and G-type for the different 
wavevectors are from Ref. 47. The Aka values agree with the 
/ 2 -JV eigenvalues in Table II of Ref. 4. 


(kx 1 ky , kz ) 

l^lx, ^y, ilz\ 

^ka 

(0,0,0) (FM, B-type) 

[100], [010], [001] 

0 

(|,0,0) (A-type) 

[100] 

-9.6874 


[010], [001] 

4.8437 


[001] - [100] 

14.5311 

(§,f,o) (c-type) 

[100], [010] 

-2.6767 


[001] 

5.3535 


[001] - [100] 

8.0302 

(§, §, |) (Neel- or G-type) 

[100], [010], [001] 

0 


three spins per unit cell. Then also taking into account 
the relation a = 2a nn , one obtains the ground-state eigen¬ 
value A = 48(2.38895) = 114.670, more than a factor of 
two larger (more stable) than the value of ss 51.3 listed 
for A( 0 ,o,o)[oio] f° r collinear FM in Table I. Thus the clas¬ 
sical MC simulations reveal a noncollinear ground state 
that is much more stable than the most stable classical 
collinear FM state. 

The results for the 2D spin lattices in Table I provide 
very useful reference points for 3D lattices, where the 2D 
results correspond to the limit c/a —> 00 . Indeed, in plots 
of Ak Q versus c/a for uniaxial 3D spin lattices below, we 
include horizontal dashed lines in the plots to observe the 
rate at which the 2D limits are approached with increas¬ 
ing c/a ratio within the calculated range 0.5 < c/a < 3. 


IV. EIGENVALUES AND EIGENVECTORS 
FOR THREE-DIMENSIONAL SPIN LATTICES 

A. Cubic Spin Lattices 

The eigenvalues and eigenvectors of the dipolar inter¬ 
action tensor for the cubic Bravais lattices are well-known 
but are presented here in modern notation for complete¬ 
ness and as a check on our calculation methods. Our 
parameters for sc, bcc and fee lattices are found to agree 
with previous results 4 and are listed in Tables II, III 
and IV, respectively, for various values of k along with 
the common magnetic structure designations. 47 Belobrov 
et al. carried out an exact calculation of the ground state 
spin configuration and energy of the sc lattice and found 
degenerate noncollinear and noncoplanar AFM ground 


TABLE III: Body-centered cubic spin lattice. Symbol defini¬ 
tions are the same as in Table II. The Aka values agree with 
the eigenvalues in Table IV of Ref. 4. 


( kx 5 ky , kz) 

[yx, [iy, /Tz] 

-^ka 

(0,0,0) (FM) 

[100], [010], [001] 

0 

N>|M 

O 

O 

[100] 

-9.6874 


[010], [001] 

4.8437 


[001] - [100] 

14.5311 

(U>0) 

[001] 

5.3534 


[110] 

7.9437 


[110] 

-13.2971 


[001] - [110] 

18.6505 


[001] - [110] 

-2.5903 

(U4) 

[100], [010], [001] 

0 

(1,0,0) 

[100], [010], [001] 

0 


TABLE IV: Face-centered cubic spin lattice. Symbol defini¬ 
tions are the same as in Table II. The designations of the 
AFM type are from Ref. 49. The Aka values agree with the 
eigenvalues in Table V of Ref. 4. 

{kx 1 ky 5 kz ) 

\Pxy fJ'Vi l*>z\ 

Akct 

(0,0,0) (FM) 

[100], [010], [001] 0 

(|,0,0) (Type-IA AFM) 

[100] 

-25.679 


[010], [001] 

12.8393 


[001] - [100] 

38.518 

(i,i,0) (Type-IV AFM) 

[110] 

14.383 


[110] 

-19.736 


[001] 

5.3535 


[110] - [110] 

34.119 


[110] - [001] 

9.029 

(0,0,1) (Type-I AFM) 

[100], [010] 

8.6687 


[001] 

-17.3374 


[100] - [001] 

26.0061 

(i,A,A) (Type-II AFM) 

[111] 

-28.9204 


[211], [Oil] 

14.4602 


[211] - [111] 

43.381 

V 3 ’ 3 ’ 3/ 

[111] 

-30.0587 


1211], [Oil] 

15.0293 


[211] - ]111] 

45.0881 

(A, 0,1) (Type-Ill AFM) 

[100] 

6.3040 


[010], [001] 

-3.1520 


[100] - [010] 

9.4560 


states with energy per spin corresponding to eigenvalue 
A = 5.344, 48 which is essentially the same as our value 
^( 1 / 2 , 1 / 2 , 0 )[ 001 ] = 5.3535 for collinear AFM ordering with 
wavevector k = (A,A,0) in Table II. The designations 
of the AFM type for fee lattices in Table IV are from 
Ref. 49. 


B. Simple Tetragonal Spin Lattices 

The eigenvalues and eigenvectors for FM moment 
alignments [k = (0,0,0)] for simple tetragonal lattices 
with c/a = 0.5-3 are shown in Fig. 21 in Appendix C and 
in a table in the Supplemental Material. 31 For c/a < 1, 
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FIG. 3: (Color online) Eigenvalues (a) Ap/ 2 , 0 , 1 / 2 ) for AFM 
wavevector k = ( 1 / 2 , 0 , 1 / 2 ) r.l.u. and (b) Aq/ 2 , 1 / 2 , 1 / 2 ) f° r 
AFM wavevector k = (1/2,1/2,1/2) r.l.u. versus the c/a ra¬ 
tio for a simple tetragonal lattice with the moments aligned 
along [ 010 ] (6 axis, solid red circles), [ 001 ] (c axis, solid green 
diamonds) and [ 100 ] (a axis, solid blue squares). 


moment alignment along the c axis is energetically favor¬ 
able, whereas for c/a > 1 , a&-plane moment alignment is 
preferred. 

The eigenvalues and eigenvectors for a number of 3D 
AFM structures for simple tetragonal spin lattices were 
determined versus c/a. The 2D limits corresponding to 
c/a —1 00 are shown as black horizontal dashed lines 
in the figures. The data for k = (1/2,0,0) are plotted 
in Fig. 22 in Appendix C. In this case there are three 
distinct Aka values for the three eigenvectors [ 100 ], [ 010 ] 
and [ 001 ] because this k breaks the fourfold rotational 
symmetry about the c axis, with the easy axis switching 
from [001] for c/a < 1 to [010] for c/a > 1. One sees 
that the respective 2D limits in Table I are reached for 
c/a > 2. Similarly, data for k = (1/2,1/2,0) and (0,0,1/2) 
are plotted in Fig. 23 in Appendix C and the data are 
listed in the Supplemental Material . 31 

The eigenvalues for AFM wavevectors k = (1/2,0,1/2) 


BaMn 2 As 2 


BaFe 2 As 2 



Ba 


Fe 

As 


FIG. 4: (Color online) Crystallographic structures of bet 
ThCr 2 Si 2 -type BaMn 2 As 2 and BaFe 2 As 2. 33 The magnetic 
atoms Mn and Fe form simple-tetragonal sublattices. The 
AFM structure of BaMn 2 As 2 is Neel-type (G-type) with AFM 
propagation vector k = (|, §, and with the ordered mo¬ 
ments aligned along the c axis, whereas the AFM structure 
of BaFe 2 As 2 is stripe-type with AFM propagation vector 
k = ( 1 , 0 , |) and with the ordered moments aligned along the 
a axis of the simple-tetragonal Fe sublattice structure (due to 
an orthorhombic distortion, the a and b axes have slightly 
different lengths at T < TW in BaFe 2 As 2 ). 


and ( 1 / 2 , 1 / 2 , 1 / 2 ) are plotted for the respective eigenvec¬ 
tors versus the c/a ratio for a simple tetragonal lattice 
in Figs. 3(a) and 3(b), respectively, with the numerical 
values listed in the Supplemental Material . 31 Here again, 
the respective 2D limits in Table I are reached rather 
quickly with increasing c/a in Fig. 3 at c/a ~ 2. 

Shown in Fig. 4 is the bet ThCr 2 Si 2 -type crystal 
structure (space group IA/mmm ) of BaMn 2 As 2 and 
BaFe 2 As 2 . 33 In both compounds the transition-metal 
atoms Mn and Fe form a simple tetragonal sublattice 
with lattice parameters ciMn/Fe = <ibct /\/2 and CM n /Fe = 
Cbct/2, yielding c Mn /Fe/aMn/F e = {c M /a hct )/s/2 = 2.285 
for BaMn 2 As 2 and 2.32 for BaFe 2 As 2 - BaMn 2 As 2 has a 
G-type (Neel-type) AFM structure with k = (i, |) in 

the tetragonal lattice notation below Xn = 625 K with 
the Mn ordered local moments aligned along the c axis, 
whereas BaFe 2 As 2 has a stripe-type itinerant AFM struc¬ 
ture with k = (A,0, 0) below Tn = 137 K with the Fe 
ordered moments aligned along the a axis of the simple- 
tetragonal sublattice in Fig. 4. The ordered moment axis 
for BaMn 2 As 2 agrees with the prediction for the wavevec¬ 
tor k = ( 3 , 5 , 3 ) in Fig. 3(b). However, as shown in 
Fig. 3(a), for BaFe 2 As 2 MDIs favor the b = [010] easy 
axis for k = ( 5 , 0 , and c/a = 2.32, perpendicular to 
the in-plane component k a b = ( 5 , 0 ) of the AFM prop¬ 
agation vector, whereas the easy axis is found to be the 
a axis, parallel to k„f, (see Fig. 40 of Ref. 33). Therefore, 
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c/a 

FIG. 5: (Color online) (a) Dependence of the eigenvalue 
A(o,o,o)[o,o,i] on the e/a ratio for a bet lattice with a FM align¬ 
ment of the magnetic moments along the c axis, (b) Expanded 
plot of the data in (a) for 0.95 < c/a < 1.5. One sees that 
FM alignment along the c axis is the most stable for c/a < 1, 
which from Fig. 21 is also the case for the simple tetragonal 
lattice. For 1 < c/a <1.3 the easy axis for FM alignment is 
the a or b axis, for 1.3 < c/a < %/2 the c axis is favored, then 
for c/a > \/2 the a or b axis is again favored. For this magnetic 
structure, one has A^o^jioo] = Afo^oiioio] — — A(o.o,o)[ooi]/2. 

there must be another source of anisotropy in BaFe 2 As 2 
that overcomes that due to MDIs to determine the easy 
axis. 


C. Body-Centered Tetragonal Spin Lattices 

The behavior of the eigenvalue A( 0 ,o.o)[ooi) of the MDI 
tensor for FM ordering with k = (0,0,0) and the ordered 
moment direction along the c axis is shown versus c/a 
in Fig. 5(a). A list of the numerical data is given in the 
Supplemental Material . 31 An expanded plot of the data 
for 0.85 < c/a < 1.5 is shown in Fig. 5(b). One sees an 
S-shaped oscillation in the latter range that was appar¬ 
ently first noticed by Lo et al. 19 The first zero crossing 



FIG. 6 : (Color online) Eigenvalues for wavevector k = 
( 5 , 5 , 0 ) r.l.u. versus the c/a ratio for a bet spin lattice with 
the moments aligned along (a) [ 0 , 1, 0 ] (solid red circles), 
[— \pY — x 2 , 0, x] (solid green diamonds) or [: x , 0, %/l — x 2 ] 
(solid blue squares), where x versus c/a is shown in (b). 

occurs at c/a = 1 , corresponding to a bcc lattice, and 
the third zero crossing occurs at c/a = \[2. This latter 
c/a value for the bet lattice corresponds to an fee lattice 
within the bet lattice that is rotated by 45° with respect 
to the bet lattice as shown in Fig. 15 of Ref. 33. The 
lattice parameters are related by af cc = y/ 2 abct = Cbctj 
yielding Cbct/abct = \/2. Hence both values c/a = 0 
and \/2 correspond to cubic Bravais lattices, for which 
it is well known that Amo.o) = 0 for all /t. We verified 
that our A(o,o,o)[ooil versus c/a data in Fig. 5(b) calcu¬ 
lated by direct summation quantitatively agree with the 
corresponding eigenvalue data in Refs. 18-20 that were 
calculated using the Ewald-Kornfeld method. 

The eigenvectors and eigenvalues of G, were calculated 
for several AFM propagation vectors. The An/ 2 ,i/ 2 ,o)[ooil 
data for k = (A, A, 0 ) are plotted versus c/a in Fig. 24 in 
Appendix C and a listing of the numerical data is given in 
the Supplemental Material . 31 The eigenvalues for wave k 
= ( ', , 0 , 5 ) versus the c/a ratio with the moments aligned 
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GdCu 2 Si 2 


4 j 



(101) plane 



FIG. 8: (Color online) Eigenvalues for wavevector k = 
(0,0,1) r.l.u. versus the c/a ratio for a bet spin lattice with the 
moments aligned along [001] (c axis, solid green diamonds) or 
[100] or [010] (a or b axis, solid blue squares). The 2D limits 
for the square lattice for c/a —> oo are shown as horizontal 
black dashed lines. 


FIG. 7: (Color online) Crystal and magnetic structures of 
bet CdCu 2 Si 2 with the ThCr 2 Si 2 -type crystal structure. One 
crystallographic unit cell is shown. The magnetic unit cell 
has dimensions 2a x b x 2c and contains four crystallographic 
unit cells. The collinear magnetic structure has an AFM 
propagation vector (1,0, 1) r.l.u. perpendicular to the (101) 
plane shown, with the magnetic moments oriented along the 
b axis. 50 Within each such (101) plane, the Gd magnetic mo¬ 
ments are FM aligned. 


in the [0, — 1 , 0 ], [— \/l — x 2 , 0, x] or [x, 0, \/l — x 2 ] direc¬ 
tions are plotted in Fig. 6(a), and x versus c/a is plotted 
in Fig. 6(b). The numerical data in Fig. 6 are listed in 
the Supplemental Material. 31 

The compound GdCu 2 Si 2 has the bet ThCr 2 Si 2 -type 
structure with space group IA/mmm as shown in Fig. 7 
and lattice parameters and z-axis Si positions a = 
3.922 A, c = 9.993 A, c/a = 2.548 and ^si = 0.368 
at 24 K. 50 The magnetic structure of GdCu 2 Si 2 is 
collinear, with the Gd ordered moments oriented along 
the tetragonal b axis with an AFM propagation vector 
k = (1,0,1) r.l.u., 50 as shown in Fig. 7. The ordered 
moment at 2 K is 7.2(4) pe/Gd, 50 in agreement with 
the value of 7/re/Gd obtained from the usual relation 
H S at = where here S = 7/2 and g = 2. Thus 

the Gd moments in (101) planes are FM aligned and are 
oriented perpendicular to k. From Fig. 6, dipolar inter¬ 
actions for k = (2, 0 , 5) and c/a = 2.548 predict that the 
moment alignment should be along the b axis, in agree¬ 
ment with the experimental AFM structure in Fig. 7. 

The eigenvalues for AFM propagation vector k = 
(0,0,1) in the bet spin lattice versus the c/a ratio with 
the moments aligned along the c axis or in the ab plane 
are plotted in Fig. 8 and listed in the Supplemental 
Material. 31 

The compound EuCu 2 Sb 2 has a primitive tetragonal 
CaBe 2 Ge 2 -type crystal structure (space group PA/nmm ) 


FIG. 9: (Color online) Crystallographic and AFM A-type 
structure with k = (0,0,1) and p = [100] of EuCu 2 Sb 2 with 
c/a = 2.401 (Refs. 29, 30) and MnF 2 with c/a = 0.6793 and 
p = [001]. j2 ~ OJ Each compound contains a bet sublattice of 
magnetic ions, but with c/a < 1 and c/a > 1, respectively, 
which is the crossover point between [001] and [100]- or [010]- 
axis ordering, respectively. 



EuCu 2 Sb2 


Mn 
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containing Eu +2 ions in crystallographically-equivalent 
sites forming a bet sublattice as shown in one panel of 
Fig. 9. 29 Like Gd +3 , the Eu +2 ions have spin S = 7/2, 
g = 2, angular momentum L = 0 and a saturation mo¬ 
ment /z sa t = gSg b = 7 /j,b- The compound orders an- 
tiferromagnetically below Xn = 5.1 K with an A-type 
structure, k = (0,0,1), and with the Eu +2 moments ori¬ 
ented in the ab plane as shown in Fig. 9. 29,30 The powder 
neutron diffraction measurements 30 can only determine 
that the ordered moments lie in the ab plane and not 
their direction within this plane. 51 EuCu 2 Sb 2 has lattice 
parameters a = 4.488 A, c = 10.778 A and c/a = 2.401. 
From Fig. 8 , for this c/a value the ordering direction for 
k = (0,0,1) is predicted for MDIs to be in the ab plane, 
in agreement with the experimental data. 

The compound MnF 2 has the primitive tetragonal ru¬ 
tile crystal structure with space group P4 2 /mnm and 
is widely considered to be the prototype for collinear 
AFM ordering. The crystal and magnetic structures of 
MnF 2 are shown in Fig. 9. At T = 298 K, the lat¬ 
tice parameters are a = 4.8734(2) and c = 3.3099(5) A 
with c/a = 0.6792 and the general F position parame¬ 
ter is it = 0.310(3). 52 The Mn +2 d 5 ions with expected 
high-spin S = 5/2 form a bet sublattice. The Mn +2 
spins order in an A-type AFM structure 53 below the 
Neel temperature In = 67 K (Ref. 54) with an or¬ 
dered moment at 5 K of 5.12(9) /ie/Mn . 55 The ordered 
moment is in good agreement with the expected value 
Msat = gSg b = 5/iB/Mn for g = 2 . A fit to x(T) mea¬ 
surements from 200 to 300 K by the Curie-Weiss law gave 
a molar Curie constant of 4.47 cm 3 K/mol and a Weiss 
temperature 0 = —97.0 K . 56 The Curie constant is close 
to the value of 4.38 cm 3 K/rnol expected for S = 5/2 and 
g = 2. From the c/a ratio and Fig. 8 , the MDI favors 
ordered moment alignment along the c axis, in agreement 
with the easy axis observed in Fig. 9. This ordering axis 
is perpendicular to the ordering axis for EuCu 2 Sb 2 with 
c/a > 1 discussed above, as expected from MDIs. 

D. Simple Hexagonal (Triangular) and Honeycomb 
Spin Lattices 

The eigenvalues and eigenvectors of the MDI tensor G,; 
for stacked simple hexagonal lattices were calculated ver¬ 
sus c/a from 0.5 to 3 for FM alignment (k = 0) and AFM 
wavevectors k = ( 1 , 0 , 0 ), (^, 0 ), (|, 3 , tt), ( 3 , 0 ) and 

(|, |), and are plotted in Figs. 25, 26 and 27 in Ap¬ 

pendix C and the numerical data are listed in the Sup¬ 
plemental Material . 31 In contrast to the AFM cases, for 
the FM alignment the approach of the eignevalues to the 
asymptotic 2D ones with increasing c/a is very slow as 
seen from comparison of the plots for FM alignments in 
Fig. 25(a) with the AFM ones, which reach their 2D val¬ 
ues by c/a ~ 2 . 

The eigenvalues and eigenvectors of G, for the hon¬ 
eycomb spin lattice in Fig. 1 calculated versus c/a from 
0.5 to 3 for k = (0,0,0) (FM alignment) and AFM prop¬ 


agation vectors k = (T, 0, 0), (0,0,^) (Neel-type in all 
directions), (0,0,0) r.l.u. (Neel-type in ab plane and FM 
alignment along c axis) and (0,0, ^) (FM alignment in¬ 
traplane and AFM alignment interplane) are plotted in 
Figs. 28, 29 and 30 in Appendix C, respectively, and are 
listed in the Supplemental Material . 31 Similar to the be¬ 
havior of the eigenvalues for the simple hexagonal spin 
lattice, for FM alignment in the honeycomb lattice the 
approach of the eigenvalues to their 2D limits with in¬ 
creasing c/a is very slow compared to behaviors for the 
AFM moment alignments. For the Neel AFM alignments 
both just in the ab plane and also along the c axis, the 
approach with increasing c/a to the infinite c/a limits is 
very fast, being essentially complete by c/a ~ 1.5. 

V. EIGENVALUES AND EIGENVECTORS FOR 
NONCOLLINEAR ANTIFERROMAGNETS 

The relationship between the ordered/induced central 
moment /A and another moment jlj at position Vji with 
respect to Jli in a collinear magnetic structure was given 
in Eq. (14). In noncollinear AFMs one must specify the 
directions of each of the moments in a crystal in order 
to calculate the net dipolar interaction of a given central 
moment /q with its neighbors inside the Lorentz sphere. 
There are two generic cases. In the first, one can define a 
nonzero AFM propagation vector k such that moments 
in a plane perpendicular to k are FM-aligned and all 
change their directions from plane to plane along k. In 
the second, the spin lattice is a non-Bravais lattice and 
the magnetic and chemical unit cells are the same, where 
the AFM propagation vector is k = (0,0,0) for such cases. 
We consider the first type of AFM ordering in the 2D tri¬ 
angular lattice in the following section and then the sec¬ 
ond type of ordering in GdB 4 and the Shastry-Sutherland 
lattice. 


A. 2D Triangular Lattice Antiferromagnets 

It is well known that the classical ground state of a 
triangular lattice AFM interacting by isotropic Heisen¬ 
berg exchange is the coplanar noncollinear 120 ° struc¬ 
ture, where each of the six neigbors of a given moment 
is at a 120 ° angle with the given moment, as in the cy¬ 
cloidal AFM structure shown in Fig. 10 where the 2D 
AFM propagation vector is k = (|,|) r.l.u. In the 
absence of anisotropy, the energy of the spin lattice in 
Fig. 10 is invariant on rotating each spin by the same 
angle, thus retaining the 120 ° angles between adjacent 
moments. Here we examine whether the MDI can de¬ 
termine how the moments are oriented with respect to 
the hexagonal unit cell axes for the AFM structure in 
Fig. 10, or indeed whether the MDI alone can stabilize 
this magnetic structure. 

The approach we use is to first calculate the eigen¬ 
values of the MDI tensor G^ for noncollinear moments 
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120° Cycloidal Ordering 


A 

k 



FIG. 10: (Color online) Coplanar noncollinear magnetic unit 
cell of classical 120° ordering on the 2D simple hexagonal 
(triangular) spin lattice for cycloidal AFM ordering with a 
commensurate wavelength of 3a/2. The hexagonal lattice 
translation vectors a and b (a = b) and the direction k of 
the cycloid wavevector k are indicated. The long-dashed line 
is the outline of the hexagonal unit cell containing one spin 
and the solid line is the outline of the magnetic unit cell con¬ 
taining nine spins (nine unit cells). The AFM propagation 
vector is k = (|,A) r.l.u. The quantity d is the distance 
between lines of FM-aligned magnetic moments along the cy¬ 
cloid axis (k) direction. The rotation angle of the magnetic 
moments between adjacent lattice lines in the k direction is 
(f>ji = kd = ^ rad. 

and variable k = ( x , x) r.l.u. and see whether the max¬ 
imum eigenvalue is obtained for x = 1/3. If so, then 
we are done. If not, we conclude that exchange interac¬ 
tions alone determine k = (|,|) r.l.u. and then calculate 
for this k what the moment orientations should be with 
respect to the crystal axes as predicted by the MDI. 

The MDI tensor Gj was calculated using Eq. (19). 
Shown in Fig. 11 are plots of the two eigenvalues \( x ,x) 
versus x with k = (x, x ) r.l.u. for the two eigenvectors 
/tii and /ti 2 shown in the figure for the orientation of 
central moment /q at the origin of the Cartesian coor¬ 
dinate system (the third eigenvalue is for FM ordering 
along the c axis as discussed in Sec. IIF and is not rele¬ 
vant here). From Fig. 11, there is no maximum in \( XtX ) 
at x = 1/3 corresponding to the 120° noncollinear struc¬ 
ture. Instead, the MDI favors k = (1/2, 1/2) r.l.u. Set¬ 
ting x = 1/3, we obtain Ao^i/s) = —1.1659 for the two 
degenerate eigenvectors /q = [100] or [010]. The AFM 
structure in Fig. 11 corresponds to Zb = [010]. 

Interestingly, the eigenvalue = —1.1659 is 

identical to the value in Table I obtained for collinear 
AM AFM ordering on the triangular lattice with k = 
A, A) r.l.u. with the same two eigenvectors. This shows 
that the net energy of interaction of a moment with the 
magnetic fields of the other moments inside the Lorentz 
sphere only depends on the projections of those moments 
on the eigenvector axis. 

The fact that is negative, whereas the eigen¬ 

value for collinear AM ordering along the easy c axis for 
k = ( 5 , 5 ) in Table I is positive, suggests that the MDI 



FIG. 11: (Color online) Variation in the eigenvalues 
versus x in the AFM propagation vector k = ( x,x ) r.l.u. for 
the two eigenvectors fin and fli 2 of the MDI tensor for the ori¬ 
entation of a representative moment fli. The first eigenvector 
is in the hexagonal b direction and the second is in the b* 
direction, which is rotated clockwise by 90° from the first (see 
Fig. 16 in Appendix A). The two curves cross at x = 1/3 and 
x = 2/3. For x = 1/3 the eigenvectors are calculated as [100] 
(along the a axis) and [010] (perpendicular to the a axis). 

might tend to cant the moments in the classical 120 ° 
coplanar structure out of the ab plane and also introduce 
an amplitude modulation of the ordered moments. 

B. GdBi and Shastry-Sutherland Antiferromagnets 

The AFM structure for GdB 4 shown in Fig. 12 
was deduced from neutron diffraction measurements . 57 
The configuration of the exchange interactions J± and 
J 2 shown in the figure is an example of a so-called 
Shastry-Sutherland Heisenberg exchange model in two 
dimensions . 58 In GdB 4 , this AFM structure is stacked 
along the c axis with FM alignments between nearest- 
neighbor layers and a corresponding FM interlayer inter¬ 
action J c that is not included in the Shastry-Sutherland 
model. 

Here we assume that the AFM structure is known, 
along with the relative orientations of each of the ordered 
moments in a unit cell. For noncollinear AFMs, Eq. (14) 
cannot be used and instead one must express each /t*, in a 
magnetic = crystallographic unit cell in terms of the cen¬ 
tral moment fli around which the dipolar sum within the 
Lorentz sphere is calculated. Thus we use the method 
described in Sec. HE to obtain the orientation (eigen¬ 
vector) of fli with respect to the Cartesian coordinate 
system, together with the associated eigenvalue. 

GdB 4 has a primitive-tetragonal crystal struc¬ 
ture with space group PA/mbm. 57 The Gd atoms 
occupy the Wyckoff 4 g positions ( 1 ) (A — x, x, 0 ), 
(2) (l — x, \ — 1,0), (3) (|+x, 1 — x,0) and 
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FIG. 12: (Color online) Four crystallographic and magnetic 
unit cells of the Gd sublattice of the tetragonal GdB4 com¬ 
pound in the ab plane. 40 The Gd ordered moments all lie in 
the afr-plane in [110] and equivalent directions. 57 Also shown 
are the 2D in-plane Shastry-Sutherland 58 exchange interac¬ 
tions Ji and J 2 between nearest- and next-nearest-neighbor 
Gd spins, respectively. The four Gd spins in the lower-left 
unit cell are numbered counterclockwise as shown. The spin 
interactions are topologically the same as in the undistorted 
Shastry-Sutherland square lattice model in which the GdB 4 
squares are not tilted with respect to the a and b axes. Ad¬ 
jacent stacked layers along the c axis are FM-aligned with 
FM (negative) nearest-neighbor exchange interaction J c (not 
shown). Since the chemical and magnetic unit cell are the 
same, the AFM propagation vector is k = (0,0,0). 


(4) (cc, \ + x, 0) with x = 0.31746(2). Thus from 
Eq. (Ala) the absolute positions of the atoms within the 
unit cell normalized to the a-axis lattice parameter are 


r i / 1 c\ 

— = [n a + - - x, n b + x, n c — , 
a \ 2 a J 


(30a) 


£2 

a 



+ 1 - x, n b + - - x, n c 



(30b) 


£3 

a 



rib + 1 — x, 



(30c) 


£4 

a 


( 1 

I n a +x, n b + -+x, 



(30d) 


where n a , n b and n c are positive or negative integers or 
zero. Taking the central moment /q to be at position ri 
with n a = n b = n c = 0, one obtains the r ^ = r/. — rj as 


£l£ 

a 

£2i 

a 

£3i 

a 

£ 4 » 
a 


( 

( 

( 


Tlai Tl b , Tl c J 1 

aJ 

1 1 0 
n a + n b + - - 2 x, 

n a + 2x, n b + 1 — 2 x, 

1 1 

n a ~ - + 2 x, n b + -, 



(31a) 

(31b) 

(31c) 

(31d) 


The 3x3 rotation matrices for the four numbered 


TABLE V: GdB 4 and Shastry-Sutherland lattice. Eigen¬ 
vectors [n x , Hy, fi z ] for central spin fli and eigenvalues Aka 
for the 3D Gd sublattice in GdB 4 and for the 2D Shastry- 
Sutherland model. For both cases, the experimental 90° 
angles between adjacent spins was assumed with the order 
(j> = <t> 0 , 4 >0 + 90°, 0o + 100° and 0o + 270° on going clock¬ 
wise around a Gd square as shown in Fig. 12, but with the 
value of 0o undetermined for the moment in the lower left 
corner of each square. The experimental x value and c/a ra¬ 
tio for GdB 4 are 0.317 46 and 0.567 97, respectively. In the 
2D Shastry-Sutherland model, x = 1/4 and c/a = 00 . The 


symbol 1 means 

-1. 



System 

X 

[fix, fly, flz] 

^ka: 

3D GdB 4 

0.31746 

[001] 

27.945 


(actual) 

[110] 

20.112 



[110] 

-48.055 



[ 001 ] - [110] 

7.883 



[001] - [110] 

76.000 


1/4 

[110] 

40.013 



[001] 

26.833 



[110] 

-66.845 



[110] - [001] 

13.180 



[110] - [110] 

106.858 

2D Shastry- 

1/4 

[110] 

40.790 982 

Sutherland 


[001] 

7.483 697 



[110] 

-48.274 678 



[110] - [001] 

33.307 285 



[110] - [110] 

89.065 660 


moments in the lower-left unit cell in Fig. 12 are 


Ri = 1) 

(32a) 

R 2 = yx-xy, 

(32b) 

R 3 = — 1, 

(32c) 


R 4 = xy-yx, (32d) 

where xy and yx are 3x3 dyadics. 

The sums in Eq. (17a) were calculated out to a Lorentz 
sphere radius R/a = 50 for 3D GdB 4 . Then diagonaliz¬ 
ing Gi gave the eigenvectors and corresponding eigenval¬ 
ues listed in Table V. Recalling that the largest positive 
eigenvalue corresponds to the minimum energy according 
to Eq. (16a), the data in Table V show that the MDI fa¬ 
vors moment alignment along the c axis, contrary to the 
experimental result in Fig. 12 which gives the alignment 
of the k = 0 spin as the [1,1,0] direction, correspond¬ 
ing to the second-highest Ak Q . The highly unstable [110] 
direction for central moment #1 corresponds to all mag¬ 
netic moments in Fig. 12 rotating clockwise by 90° and 
hence all moments in each Gd 4 square pointing towards 
the center of the square. The RKKY interaction between 
Gd spins and/or a high-order crystalline electric field ef¬ 
fect evidently give an anisotropic exchange interaction 
that is responsible for the observed ordered moment di¬ 
rections. 

Calculations were also carried out for x = 1/4, which 
corresponds to untilted Gd 4 squares in Fig. 12, as shown 
in Table V. One sees significant differences in the 
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eigenvalues compared to the results for the observed 
x = 0.31746. In particular, the Gd ordered moments 
are now predicted to have the experimental ordered- 
moment directions. We also carried out calculations for 
the Shastry-Sutherland 2D lattice and the results are 
shown in Table V, where the favored ordered moment di¬ 
rection for Gdi is found to be the same as for x = 1/4 and 
c/a = 0.56797, the observed c/a ratio for GdB 4 . Thus 
the ground-state ordering direction predicted by the MDI 
is sensitive to the tilting angle of the Gd 4 squares. 


VI. MAGNETIC ORDERING TEMPERATURE 
DUE TO MAGNETIC DIPOLE INTERACTIONS 

The MFT calculations in this and the following sec¬ 
tions closely follow the development of the author de¬ 
tailed in Ref. 11. Therefore only an outline of the calcu¬ 
lations associated with the MDI is given. 

In this section, an AFM ordering (Neel) temperature 
arising from dipolar interactions only is denoted by Tna 
and a FM ordering (Curie) temperature by Tca, where 
the subscript A refers to the quantity being the contribu¬ 
tion from an anisotropic magnetic interaction. Similarly, 
a Neel temperature arising from Heisenberg exchange in¬ 
teractions only is denoted by Tnj and a Curie tempera¬ 
ture by Tqj. We use the Weiss MFT to calculate these 
transition temperatures where we assume that the spins 
are identical and crystallographically equivalent and we 
only treat EM (not AM) magnetic structures on Bravais 
lattices. Within MFT, the contributions of the dipolar 
and exchange interactions to the actual ordering temper¬ 
atures Tn and Tc, respectively, are additive: 

Tk = T na + T N j, Tc=Xca+7cj- (33) 

The magnetic ordering temperature T m j (m = N, C) 
for both AFMs and FMs due to exchange interactions is 
given by the same expression 11 


The magnetic ordering temperature is determined 
within MFT by the criterion that /2 0 —* 0 for T —> T~. 
For magnetic dipole ordering, the near-field contribution 
to the local magnetic induction is given by Eq. (20). The 
magnetic moment ft in that equation is defined in gen¬ 
eral as either the ordered moment in a magnetic struc¬ 
ture in H = 0 (/io) and/or an average moment induced 
by H a > 0 (/x). Using Eq. (35), Eq. (20) associated with 
MDIs becomes 


td near 
-^int o: i 


9 MB S M-^ka: 


(37) 


A. Antiferromagnetic Ordering (Neel) 
Temperature 

Here we calculate Tna in H = 0 within MFT for a 
specified AFM wavevector k and ordered moment axis p 
in the presence of MDIs but in the absence of exchange 
interactions. The standard MFT prediction is obtained 
from 11,40 ’ 59 


go = Bs 


V k B T )' 


(38) 


where we have dropped the subscript i because all mo¬ 
ments are crystallographically equivalent in H = 0, the 
subscript 0 in po signifies that H = 0 as above, and 
Bs(y) is the Brillouin function for spin 5 given by our 
unconventional expression 


Bs(y) = 


^{< 2S + 1 > 


coth 


(25+1) 


— coth 


(I)}- 


(39) 

There is no demagnetizing field for an AFM in H = 0 be¬ 
cause there is no net magnetization, so for AFM ordering 
in H = 0 the local field is just the near field. Inserting 
T^inTai from Eq. (37) into (38) gives 


T m j = — 


5(5+1) 

3/cb 


Y Jij cos 


(34) 


where (f>ji is the angle between magnetic moments j and i 
in the ordered state and (j)ji = fa — fa = 0 for a FM. 
We define the reduced ordered and/or applied magnetic 
field-induced average moment p for a spin 5 as 


jj_ = ft 

ftsat Z?5fi B 


(35) 


where ft sat = gSg b is the saturation moment of the spin 
and g is the spectroscopic splitting factor. Using Eq. (34), 
one can write the exchange field seen by a representative 
moment i in zero applied field H as 


ZJ _ Tmj 

H exchi — MO 

W 


3k B T mJ _ 
5Mb(5 + 1) M °’ 


(36) 


where the subscript 0 in p 0 signifies H = 0, C\ is the 
single-spin Curie constant (see below) and this expression 
applies to the ordered state. 


where 


fto = B s (yo), 


V o 


g 2 Sg 1/ipAka 
a 3 k B T 


(40a) 


(40b) 


Then one obtains for a given k and easy axis a the Neel 
temperature 11 


Tnao 


g 2 5(5 + l)ft|A ka 
3a 3 k B 


(41) 


The relevant ordering axis a and hence Tnaq is the one 
with the largest eigenvalue Aka for the given AFM struc¬ 
ture. 

The single-spin Curie constant C\ for spin 5 is given 
by 37 


, = g 2 5(5 + l)/4 
3 /cb 
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so Eq. (41) can be written more succinctly as 


TNAa 


Ci A 


l^k a 


Thus one can also write 

^kc* _ ^NAa _ 3fcB^NAa 
a 3 Cl g 2 S(S + l)/x|' 


(43a) 


(43b) 


Then for H = 0 and T < Tna one can write the near 
field in Eq. (37) in the direction of each ordered moment 
as 


Rnear _ 3fc B T NA a/ro 
inta g(S + l)g B ' 


(44a) 


The exchange field for H = 0 seen by each moment in 
its ordering direction due to Heisenberg exchange inter¬ 
actions for either FM or AFM ordering can be written in 
the same form as 11 


r^exch 


3A'ii / m .//Ao 

g{S + 1)/a b 5 


(44b) 


where T m j is the contribution of Heisenberg exchange 
interactions to either a FM Curie temperature Tqj or 
an AFM Neel temperature Tnj. Using Eq. (33), in the 
case of AFM ordering the sum of the two local fields in 
Eqs. (44) can be written 


nlocal _ 3fc B (r N J + T N Aa)Ah) _ 3fc B T N /rQ 
g(S+1 )hb g(S+ l)/i B ’ 


(45) 


where Tn is the Neel temperature in the presence of both 
exchange and MDIs. 

Because different sources of local fields are additive in 
their contributions to the observed Tn within MFT, if 
both exchange and dipolar interactions are present Tna 
is the contribution of dipolar interactions to Tn, which is 
usually but not always a small fraction of Tn- 

Quantum fluctuations generally increase as the dimen¬ 
sionality of a spin lattice decreases. These quantum fluc¬ 
tuations can prevent long-range magnetic ordering from 
occurring. Corruccini and White found from spin-wave 
calculations that AFM order cannot occur at finite tem¬ 
perature on the 2D square spin lattice due to dipolar 
interactions alone. 13 MFT does not take into account 
such quantum fluctuations associated with reduced di¬ 
mensionality and hence predicts that AFM ordering can 
occur in ID, 2D and 3D spin lattices. 


where a is the a-axis lattice parameter of the unit cell, 
Ugpin is the volume per spin, Aoa refers to the FM mo¬ 
ment alignment, the magnetic moment per unit volume is 
t</Uspin = ggsSjl/Vs pin, and we used Eqs. (35) and (37). 
Then following the same development as in the previous 
section gives the Curie temperature 


Tca« 



(FM), (47) 


where C\ was defined in Eq. (42). The system will choose 
the easy axis a with the largest value of Aoa- For a cubic 
Bravais lattice Ao a = 0, so there is no preferred easy axis 
for FM ordering according to the present treatment. 

Using Eq. (47) one can write the local field in Eq. (46) 
for FM moment alignments as 


olocal _ 


3fc B TcAet _ 
gg B (S + 1) 


(48) 


If Heisenberg exchange interactions are present, one adds 
the local exchange field in Eq. (44b) to the dipolar con¬ 
tribution in Eq. (48) to obtain 


olocal 3fc B T Ca 
intai ~ gg B (S + 1) M 


(49) 


where Tq u = Tq A a + Tcja according to Eq. (33). 

Comparing Eqs. (45) and (49) one sees that the same 
form of the local field in the direction of each ordered 
moment is obtained for both FM and AFM structures 
in the ordered states and one can therefore write the 
local magnetic induction seen by each moment in general 
for either FM or AFM moment alignments and dipolar 
and/or Heisenberg interactions as 


olocal _ 


3 fc B T, 


me* _ 

/c 


gg B (S + 1 ) 


(50) 


where T mQ is the Curie or Neel temperature for the 
collinear ordering axis a. 


VII. COMPETITION BETWEEN 
FERROMAGNETIC AND 
ANTIFERROMAGNETIC ORDERING 


B. Ferromagnetic Ordering (Curie) Temperature 


As is well-known, whether or not a particular sample 
exhibits FM ordering driven by the MDI depends on the 
shape of the sample via the demagnetizing field as well as 
the competition with AFM states. The former is evident 
from Eq. (8a) which for H a = 0 becomes 


olocal gI 1 1' \ 

■^int ai ^3 / 'Oa: 


4t r /I 

U spin /a 3 \3 


Ada 


(46) 


One can have a crossover between FM and AFM or¬ 
dering depending on the value of the demagnetizing fac¬ 
tor Nda and the possible AFM eigenvalues Aka and FM 
eigenvalues Ao a ■ The value of Ada depends on the shape 
of the sample. For FM ordering, the field direction with 
the smallest value of Ad« gives the lowest free energy and 
hence is the FM ordering direction provided that the cal¬ 
culated TcAa > 0 and that competing AFM states have 
a lower calculated Tnao > 0 . 
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To examine this competition, we define the dimension¬ 
less reduced AFM and FM ordering temperatures ob¬ 
tained from Eqs. (43a) and (47), respectively, as 


7 N A a C 

c 1 

TcAaCt 3 

c\ 


Aka (AFM), 

. 4t r /I AT \ 

Aoa + T/ / 3 ( o da 1 


(51a) 

(FM). 

(51b) 


As an example, we consider the competition between 
FM and AFM ordering due to dipolar interactions on sc, 
bcc and fee Bravais lattices, which have Aoa = 0 and 
Vspin/a? = 1, 1/2 and 1/4, respectively. The reduced 
Curie temperature in Eq. (51b) is plotted versus N c \ a for 
sc, bcc and fee Bravais spin lattices in Figs. 13(a), 13(b) 
and 13(c), respectively. Using Eq. (51a) and the data in 
Tables II, III and IV, AFM Aka values are plotted for the 
most stable (positive) Aka value for each k as horizontal 
lines for the sc, bcc and fee lattices in Figs. 13(a), 13(b) 
and 13(c), respectively. One sees from Fig. 13 that for 
the magnetic structures considered, the ground state of 
the sc lattice is AFM-ordered with k = ( 2 , 5 , 0 ) r.l.u. 
and ordering axis jl = [ 001 ] for all values of Ada, the 
bcc lattice is unstable to FM ordering only for N da « 0 
and the fee lattice for Ada 1$ 0.03. These inferences are 
consistent with early results . 4 A sample with the shape 
of a long thin needle with the magnetization directed 
along the axis of the needle has a demagnetizing factor 
A d a « 0. 


VIII. PROPERTIES OF THE 
MAGNETICALLY-ORDERED STATE 

A. Ordered Moment and Magnetic Heat Capacity 

For either an AFM or FM with Heisenberg and/or 
MDIs, Eq. (50) gives the same form of the local magnetic 
induction seen by each spin in its ordering direction for 
T < T m . Using Eq. (50), the behavior of /2 versus t is the 
same as for pure Heisenberg interactions and is shown for 
several values of the spin S in Fig. 10 of Ref. 59. 

The magnetic energy per spin is given by 

£ma gi = (52a) 

where the factor of 1/2 derives from the fact that B] ocal 
is attributed to the neighbors of fit whereas the energy 
is equally shared by pairs of interacting spins. Inserting 
B] ocal f rom Eq. (50) into (52a) for a mole of spins with 
A = Aa where Aa is Avogadro’s number, one obtains 

r _ 3RS r]1 _ 2 i'cov'i 

77 ma g — 2(S+l)^ nl ^ 0, (52b) 

where R = Aa&b is the molar gas constant. Then the 
magnetic heat capacity C mag per mole of spins is obtained 



FIG. 13: (Color online) Reduced magnetic ordering temper¬ 
ature TmagU 3 /C\ versus the demagnetizing factor Ada with 
0 < Ada < 1 for pure magnetic dipolar ordering in (a) sc, (b) 
bcc and (c) fee Bravais spin lattices as predicted for FM and 
AFM ordering by MFT via Eqs. (51). An ordering wavevector 
is labeled as (mi, m 2 , m 3 ) r.l.u. and the ordered moment axis 
as [fj, x , fly, ii z \ in Cartesian coordinates. Values of T mag < 0 
are unphysical. For the fee lattice, the most stable AFM 
wavevector shown is k = (1/3,1/3,1/3) r.l.u. 
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as 11 - 59 

C'mag _ 3SfJ,o(t) d.Ho(t) , x 

R 5 + 1 dt~’ 1 ° ' 

where t = T/T mag and the reduced ordered moment ver¬ 
sus temperature /2 q (i) in H = 0 is obtained as described 
in Ref. 11. This equation is identical to that obtained for 
pure Heisenberg interactions, where plots of C mag /R ver¬ 
sus t for several values of S are shown in Fig. 11 of Ref. 59. 
For quantum spins, C mag decreases exponentially to zero 
for t —> 0, whereas for classical spins C mag /R -+ 1 for 
t -A 0. 


B. Dipolar Anisotropy of Uniaxial 
Antiferromagnets in the Ordered State 

Here we calculate the dipolar anisotropy of the free 
energy between EM orthogonal principal collinear mag¬ 
netic ordering axes denoted as the a and /3 axes. We con¬ 
sider collinear AFMs with noncubic spin lattices contain¬ 
ing identical crystallographically-equivalent spins. The 
lowest-order expression for the anisotropy free energy per 
spin Fj is given by the usual expression 

Fi = K\ sin 2 9, (53) 

where 9 is defined as the angle between the ordered mo¬ 
ment axis and the a axis. We derive an expression for 
K\ associated with the anisotropic MDI in terms of the 
eigenvalues and eigenvectors of the MDI tensor. 

The orientation of a representative T-dependent or¬ 
dered moment pi in the a-/3 plane in I? = 0 with 
Mo = \Pi\ is 


pi = /io(cos0d + sin 0/3), (54) 

where /iq is the T-dependent ordered moment in H = 0 
and 9 = 0 corresponds to pi parallel to the a axis. The 
corresponding T-dependent internal local field is 

B*nt cos 0 « + sin 9 (55) 

where the expression for B{°£^A is given in Eq. (8a) with 
H a = 0. The differential dFi of the magnetic free energy 
of the moment is 


dF i = ~p i -dB^f, (56) 

where the factor of 1/2 is present because arises 

from the neighboring moments of pi whereas the free 
energy per moment is equally shared between each pair 
of moments. Inserting Eqs. (54) and (55) into (56) gives 

dFi = - B\° n $i) sin(9 cos 9d9. (57) 

Integrating dFi from 9 = 0 to 9 yields 

Fi = ^{Btlti - B^fi) sin 2 9. (58) 


This expression for Fi applies to moments along the 
collinear ordering axis with angles of either ±0 to the 
a axis because the sine function is squared. Comparing 
Eq. (58) with (53) gives the anisotropy parameter K\ as 

R i = ( 59 ) 

For an AFM in the ordered state, one has B}°^C = 
Bf/fU. Inserting Bf/fU in Eq. (20) into (59) gives 

K x = £|(A kQ -A k/3 ). (60) 

From Eqs. (53) and (60) one obtains 

p. f° ( 0 = °) 

* \^?(A kQ -A k/3 ) (0 = tt/2) ' 

Therefore if A ka — A k ^ > 0, the minimum free energy 
occurs if the moments are aligned along the a axis (9 = 0) 
and hence the easy axis is the a axis, whereas if A ka — 
A k /3 < 0, the p axis (9 = 7r/2) is favored for the ordering 
axis over the a axis. These results are consistent with 
expectation because one expects a moment pi to line up 
along the axis with the largest value of in Eq. (20), 
i.e., with largest value of A k . 


C. Perpendicular Magnetic Susceptibility of 

Collinear Antiferromagnets in the Ordered State 

The Heisenberg exchange Hamiltonian has no intrinsic 
magnetic anisotropy to determine the directions of the 
ordered moments in the ordered state with respect to 
the spin-lattice axes. In this paper the only source of 
magnetic anisotropy is the MDI, and in this section we 
only consider collinear magnetic ordering. The easy axis 
is the eigenvector of the interaction tensor Gj(k) that 
corresponds to the largest eigenvalue for the given AFM 
propagation vector. 

The single-spin magnetic susceptibility x is rigorously 
defined as x = liniir-^o n{H)/H where ^ is the thermal- 
average moment of a spin in the direction of H that is 
induced by H. Here we take the easy axis to be the x axis 
and the applied infinitesimal field to be along a z axis, 
perpendicular to the x axis. The magnitude of each or¬ 
dered moment in zero field is which is T-dependent 
as shown in Ref. 59. In the presence of the perpendicu¬ 
lar field, the magnitude of the moment does not change 
in the AFM phase 4,11 and the induced moment acquires 
a component along the z axis. Including the applied in¬ 
finitesimal perpendicular field and both the exchange and 
dipolar fields and setting the net torque on a represen¬ 
tative moment equal to zero following the procedure of 
Ref. 11 yields the perpendicular susceptibility 

= _ Ci _ 

^ (Fnj + F NAx — T C az) — Opj 


(62) 
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The T-dependent ordered moment g o canceled out, so 
X_l is independent of T for T < T N , as also obtained for 
pure Heisenberg spin interactions. 11 

Several special cases occur for Eq. (62). If exchange 
interactions are negligible, the pure magnetic dipole pre¬ 
diction is obtained by setting Tnj = 0 p j = 0, yielding 


C i 

TnAx ~ TqA z 

a 3 

Akz - A 0l - 3V **/ a s 


(63a) 

(63b) 


For cubic Bravais spin lattices for which Ao a = 0 for all a , 
Eq. (63b) gives 


AT = i-■ (64) 

Akx 3V^ pin / a 3 

This result agrees, e.g., with \± obtained from the equa¬ 
tion between Eqs. (29) and (30) in Ref. 4 which includes 
in the denominator of Eq. (64) the ground state eigen¬ 
value A ka: = A(!/2,i/2,o)[ooi] = 5.351 (/ 5 in their notation) 
for the sc dipolar AFM, in good agreement with our value 
of 5.3535 in Table II. 

When dipolar interactions are negligible, Eq. (62) gives 
for the pure Heisenberg exchange model 

Xx = T C \ ( T < T N j), (65a) 

Jnj- tipj 

in agreement with Ref. 11. In the PM state at T > Tnj, 
the isotropic susceptibility per spin is given by the Curie- 
Weiss law 11 


x = ( T ^ T nj)- ( 65b ) 

Comparing Eqs. (65a) and (65b) gives 

XX = x(2nj) (T < T nj ). (65c) 


D. Perpendicular Critical Field 

As the perpendicular field is increased from zero at 
T < Tn, the induced perpendicular moment g± increases 
as 


Thus one obtains 


ffc(T) = Mo (T) 
H c (0) Mo(0) 


Mo CO 

1^ sat 


Mo (T), 


( 68 ) 


where is plotted versus t = T/Tn in Ref. 59. Since 
within MFT go{T) depends on the spin S of the moment, 
so does Near t = 1, one obtains 


tir “ v ' rr * ( ‘- 1_ > - (69) 

Previous classical calculations (not utilizing the Weiss 
MFT and hence not the Brillouin function for quantum 
spins) yielded the behavior in Eq. (69) for the whole 
temperature range 0 < t < 1, with the proportional¬ 
ity replaced by an equality. 39 In that case, expanding the 
right-hand side of Eq. (69) in a Taylor series about t = 0 
gives the linear dependence = 1 — f {t <C 1) instead 

of the exponential approach to unity for t —> 0 obtained 
for quantum spins. 


IX. CURIE-WEISS LAW IN PARAMAGNETIC 
STATE 


In the PM state above the Neel or Curie temperature, 
all moments are aligned in the direction a of the magnetic 
field H a applied along a principal axis of the spin lattice 
[the magnetic propagation vector is k = (0,0,0) = 0]. 
For Heisenberg exchange interactions, the exchange field 
in the PM state is isotropic and given by 11 ’ 40 


XZexch 


3k B 0pJ _ 

gH B {S + 1 )^’ 


(70) 


where /i is the normalized moment induced by H a and 


S(S + 1) 
pJ 3 k B 




(71) 


is the contribution to the Weiss temperature in the Curie- 
Weiss law due to Heisenberg exchange interactions. Then 
adding H exc h and H a to the local dipolar field for H a = 0 
in Eq. (46) gives the total local field seen by each moment 
as 


MX = Xx H, (66) 

where x_l is given by Eq. (62). When /i± reaches the or¬ 
dered moment go (T), the induced moments become par¬ 
allel to H and the system enters the PM state in a second- 
order transition. 4-11 Setting g± = go with increasing H , 
the critical field H c at which this happens is defined by 
Mo = X±H C , yielding 

HJT) = (67) 

Xx 


R 1 


local 


= H, 


3fc B fl P j 
gg B (S+l) 
gg B S 


M 


(72) 


A () fl 


47T 


3P sp i„/a 3 


where we assume that the demagnetizing field has been 
corrected for in experimental data and hence the demag¬ 
netizing factor Nda does not appear in this expression. 
To include it, replace the multiplicative factor | in the 
last term by | — Nd a . 
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Analogous to Eq. (38) for H a = 0, in the present case 
one has 


9 = Bs 


\ k B T ) ' 


(73) 


Inserting from Eq. (72) into (73), Taylor expand¬ 

ing the Brillouin function Bs(y) to first order in y, solving 
for jl and using Eq. (35) gives the Curie-Weiss law 


C x 

/Cot. rj-f x 

9p a = 9 pJ T @pAo^^ 


(74a) 

(74b) 


where the single-spin Curie constant C\ is given in 
Eq. (42), dp j is given in Eq. (71) and the magnetic dipole 
contribution 9 pAa to the Weiss temperature is 

dpAa = ^ ( A °“ + 3 W^) ' (?4C) 

A comparison of Eq. (74c) with (47) shows that the 
contributions of dipolar interactions to the Weiss temper¬ 
ature and the Curie temperature of a FM are the same, 
i.e., 


9pAct = TqAoli (74d) 


Eq. (74a) becomes 

- (T > T NAx ), (76) 

J- UpA a 

where 9 P Aa is given by setting the second term in 
Eq. (74c) to zero for a spherical sample, yielding 


9pA 


ol 


CiAq 


(77) 


This would be zero for a cubic Bravais spin lattice be¬ 
cause in that case Ao a = 0 for all a. The Neel tempera¬ 
ture in Eq. (43a) for the easy x axis is 

Tnax = (78) 

a 6 

and we define the ratio f Aa as 


/a« 


^pAa 

TnAx 


Aoa 

Akx 


(79) 


where the subscript A in f Aa signifies that the value of / 
arises only from the anisotropic MDI and a can be any 
of the principal axes x , y or z. 

Using Eqs. (78) and (79), the Curie-Weiss law (76) for 
a single spin can be written in dimensionless form as 


which is the same result as obtained from MFT for a 
system of local moments exhibiting a FM transition and 
interacting by Heisenberg exchange only . 11 

On the other hand, a comparison of Eqs. (43a) and 
(74c) shows that in general the contribution of dipolar 
interactions to the Weiss temperature for AFMs is not 
equal to the negative of the dipolar Neel temperature in 
Eq. (43a), as is also found in general for local-moment 
Heisenberg AFMs . 11 Thus the ratio / = 9 p /Tc for a FM 
within MFT is 


/ = 1 (FM), (75a) 

whereas in general for an AFM it is 

f a = ^ = H pAa + 9pJ < 1 (AFM). (75b) 


X. ANISOTROPIC MAGNETIC 
SUSCEPTIBILITY OF A SPHERICAL SAMPLE 
OF A PURE DIPOLAR ANTIFERROMAGNET 

In the following, we assume that the sample is in the 
shape of a sphere, which cancels the Lorentz field within 
the Lorentz cavity according to Eq. ( 8 a) and hence ame¬ 
liorates the competition of FM with AFM ordering. 


Xck 7 \ \ x 


1 


7a — /aq 


(7a > 1), 


(80a) 


where the reduced temperature t A is defined as 


7a 


T 

TnAx 


(80b) 


Note that Eq. (80a) is a law of corresponding states for 
all quantum spins S, since S only appears in C\. 

The reduced PM susceptibility at Tnax from the Curie- 
Weiss law (80a) is then 


X«(7a = 1 + ) Tnax 1 


Ci 1 - f Aa 

From Eqs. (80a) and (81) one obtains 

X«(7a) _ 1 /a« 

Xct(t A = 1 + ) t A -f Aa 

which yields the identity 

X a (t A = 1 + ) = 

Xol (7 A = 1 + ) ’ 


(t A = 1 +). (81) 


(t A > 1 ), (82) 


(83) 


as required. 


A. Paramagnetic State 


B. Perpendicular Susceptibility in the 
AFM-Ordered State 


For a dipolar collinear AFM at T > T^ Ax where the 
easy axis is defined as the x axis, the Curie-Weiss law in 


In the AFM state at T < T^ A of a strictly dipolar 
AFM, one sets Tnj = 9 p j = 0 and for spherical samples 
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Eqs. (63) yield 


where 


X±(T < ^NAx) 


Akx — A 02 

Ci 

TnAx — Tcaz ’ 


(84a) 

(84b) 


where, as above, the x axis is the easy axis for the 
collinear AFM ordering, XnAx is the associated Neel tem¬ 
perature and the 2 axis is perpendicular to the x axis, i.e., 
Xx = Xz- One can write Eq. (84b) in dimensionless form 
as 


XX^NAx _ 1 

Cl 1 - r z 


(t A < 1, z ± x), 


(85a) 


where according to Eq. (43a) and Eq. (47) modified for 
a spherical sample one has 


V 0 


f/MB MoAkx 
&b T a 3 


(88b) 


Here Bg(y ) is the Brillouin function in Eq. (39) and 
B' s (y 0 ) = [dBs(y)/dy]\y =yo . Note that the parallel sus¬ 
ceptibility in the dimensionless form in Eq. (87a) still 
depends on S since the Brillouin function on the right- 
hand side does. This contrasts with the dimensionless 
forms of the Curie-Weiss and perpendicular susceptibil¬ 
ities above for dipolar interactions that do not depend 
on S. 

Useful limits are 


r*(f A -t 0) = 00 , r*(t A —t 1) = 1, 


(89) 


yielding 


r z 


TqAz 

?NAx 


Aq z 

Akx 


Using Eqs. (81) and (85a) one obtains 
X-l(^A < 1) _ 1 - /ao 

Xo{t A = 1 + ) 1 - r z 


(85b) 


( 86 ) 


X||2nAx 

Ci 

XH^NAx 


— 0 (f A —> 0) 
1 


(fA ~t 1 )■ 


C 1 1 - f Ax 

The latter X|| expression is identical with 


(90a) 

(90b) 


Comparing Eqs. (83) and (86), one sees that in general 
the hard-axis Xz is continuous on cooling below Tnaz, 
where xx = Xz below Tnax- If Ao a = 0 for all a as in 
cubic Bravais lattices, xx is obtained for all axes below 
TnAx- 


Xx TnAx _ 

Cl ' 1-fAx 


(fA = 1 + ) 


(91) 


obtained from Eq. (81) for the Curie-Weiss law at f A = 
1 + for the held applied along the x axis. Thus X|| = Xx 
for f A < 1 joins continuously with Xx f° r fA > 1- 


C. Parallel Susceptibility in the AFM-Ordered 
State 


D. Example 


When an infinitesimal Held H = Hi is applied in the 
positive x direction along the collinear AFM ordering 
easy axis at a temperature 0 < T < Tnax, an ordered 
moment initially pointing parallel (antiparallel) to H in¬ 
creases (decreases) slightly in magnitude, where the vec¬ 
torial change d/1 = dyi is the same for both moments. 
Therefore in this section we only consider the change 
in the rr-axis component of a representative moment jli 
pointing towards the positive x axis due to the applied 
Held. 

Following Ref. 11 we obtain the dimensionless equation 


X || -In Ax 


1 


Cl T *(t A ) — /Ax ’ 
where f A is defined in Eq. (80b) and 


(f a) = 


(S + l)f A 


fAx = 


6■ 


pAa? 


3 B' s (y 0 ) ’ Tnax 

/x 0 (T) is obtained by numerically solving 


^kai 


(87a) 


(87b) 


(88a) 


As an example, we consider the simple tetragonal Bra¬ 
vais spin lattice with c/a = 0.8, 1.0 and 1.2 and AFM 
propagation vector k = (A,-1,0) r.l.u. for temperatures 
both above and below the Neel temperature. Recall that 
for /aq, the x and a axes are the easy principal axis 
for AFM ordering and any of the three principal axes, 
respectively, whereas for r 2 , the z axis must be an axis 
perpendicular to the x axis. In a real material, one must 
identify x, z and a with the appropriate crystal axes. 

The eigenvalues and eigenvectors of the dipolar in¬ 
teraction tensor taken from tables in the Supplemen¬ 
tal Material 31 are shown in Table VI along with the 
respective values of / Aq and r z defined in Eqs. (79) 
and (85b). One sees that the AFM state is stable against 
the FM state below Xnax for both c/a values, but that 
the anisotropy in the PM state at T > Tna changes sign 
between the two c/a values. 

Using the data in Table VI, Eq. (86) yields x = 

1 for the easy a axis for both c/a = 0.8 and 1.2. For 
the sc lattice with c/a = 1, one has Ao a = /a = r a = 
dpaa = 0 for all a. Therefore x(^) follows a Curie law 
for t A > 1. Also, there is no restoring force for keeping 
the easy axis parallel to the field, so the magnetization 


Mo = gyBSB s {y 0 ), 
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TABLE VI: Eigenvalues A and eigenvectors /t = [fi a , gb, g c ] of the dipolar interaction tensor for simple-tetragonal spin lattices 
with c/a = 0.8 and 1.2 and wavevectors (mi, m 2 , m 3 ) r.l.u. The data were taken from tables in the Supplemental Material. 31 
The largest eigenvalue for k = 0 is labeled as Ao a . For k = (A, 1, 0) r.l.u. the maximum eigenvector is denoted as Aka; and the 
value for the perpendicular direction as Akz- For each k, the values of /a and r z are listed as defined in Eqs. (79) and (85b), 
respectively. According to Eqs. (80a) and (85a), the parameter f Aa is relevant for the PM T range and r z is relevant for the 
AFM-ordered T range. In the table, the assignments of the x and 2 Cartesian axes to the c and a crystal axes, respectively, 
are shown in the subscripts to the parameters. 


c/a 

k (r.l.u.) 

[100] 

[001] 


/a 

r z 

0.8 

(0,0,0) 

(U.o) 

Aoz.a = -1.9691 

Akz, a = -4.7977 

Aox.c — 3.9382 
A k *,c = 9.5955 

fAx,c 

= 0.4104, / A z,a = -0.2052 

r z , a = -0.2052 

1.2 

(0,0,0) 

(U’O) 

Aoz.a = 0.9364 

Akz,a = -1.8691 

Aoz.c = -1.8728 
Akx.c = 3.7381 

/Ai,c 

= -0.5010, f Az , a = 0.2505 

r z , a = 0.2505 


flops to the perpendicular orientation whenever this is 
attempted. Thus only x_l(T) = xt^NAx) is measured for 
t A < 1. 

Shown in Fig. 14 are plots of x ^ A ^ NAx versus t A for 
c/a = 0.8, 1.0 and 1.2 illustrating the progression of the 
anisotropy in x as c/a traverses the sc of unity. To our 
knowledge no theoretically-predicted behaviors such as 
in Figs. 14(a) and 14(c) have appeared before in the lit¬ 
erature. 


XI. ANISOTROPY OF MAGNETIC 
SUSCEPTIBILITY OF A HEISENBERG 
PARAMAGNET DUE TO MAGNETIC DIPOLE 
INTERACTIONS 


In this section we assume that demagnetizing fields 
have been corrected for in experimental data and hence 
the demagnetizing factor TVda does not appear. 

In the PM state above In, according to Eq. (74c) the 
anisotropy in x can only arise from a difference in the 
dipolar Weiss temperatures along different principal axis 
directions a and /3, given by Eq. (74c) as 


Ci 

^pAa ^pA/3 — jy(Aoa Ao/?). (92) 


For cubic Bravais lattices, one has no dipolar anisotropy 
in the PM state because Ao Q = 0 for all a. Here we follow 
the approach of Keffer. 60 

For two susceptibilities Xa and X/3 measured along the 
a and /3 principal axes, one has the identity 


If the dipolar anisotropy in 6 is small compared to the 
measured average Weiss temperature 9 p , one can define 
the geometric-mean susceptibility x = \/XaX@ an d use 
Eq. (92) to obtain 

y 2 

Xa - X/3 = ^-(Aoa - A 0 / 3 ). (95) 

Here the Curie-Weiss x’ s are per spin and a is the a- 
axis lattice parameter for the particular Bravais spin lat¬ 
tice considered. The susceptibility difference per mole of 
spins is obtained by multiplying each \ on the left side 
of Eq. (95) and one x on the right by Avogadro’s number 
N a and Eq. (95) yields the molar susceptibility difference 

XMa(T) - XM/3 (T) = (Aoa - A 0/3 ). (96) 

l\ A a° 

Here we apply Eq. (96) to the primitive tetragonal 
rutile-structure collinear AFM MnF 2 with Tn = 69 K, 
which is often considered a prototype for collinear AFM 
ordering. This compound contains a bet sublattice of 
Mn +2 cations with spin S = 5/2 and an expected g — 2 
and orders into a A-type AFM structure with AFM 
wavevector k = (0,0,1) as shown in Fig. 9. The lattice 
parameters are 52,61 

a = 4.8734(5) A, c = 3.3103(10) A, -= 0.6793(3). 

a 

(97) 

For the given c/a ratio and FM k = 0 we find Ao[ooi] = 
4.3219 and A 0 [ioo],[oio] = —2.1609 yielding 

Aq[ooi] — A 0 [ioo] = 6.4828, (98) 


J__1_ = Xa - X/3 

X/3 Xa XaX/3 


or 


Xa - X/3 = XaX/3 




Using Eqs. (74), Eq. (93b) yields 


Xa - X/3 = 


Xa.X/3 

Cl 


(0 P 


0p/3)- 


(93a) 


(93b) 


(94) 


whereas for the ordering wavevector k = (0,0,1) r.l.u. 
we obtain A(ooi)[ooi] = 13.8639 and A( 0 oi)[ 100 ],[oio] = 
—6.9319, with A( 0 oi)[ooi] — A(ooi)[ioo],[oio] = 20.7958. 
These values show that the [001] moment direction is 
energetically favored by the MDI both above and below 
Tn, in agreement with experiment as follows. 

The anisotropic x(T) of MnF 2 crystals is shown in 
Fig. 15(a). 62,63 Above T N , x is found to be nearly 
isotropic. Below Tn, the data are a textbook example 
of the anisotropy expected for collinear AFM ordering, 
where in this case the easy axis is the c axis. According 












28 



FIG. 14: (Color online) Anisotropy of the magnetic suscep¬ 
tibilities \ a and Xc due to MDIs versus reduced tempera¬ 
ture tA = T/Tnax for a simple tetragonal spin lattice with 
(a) c/a = 0.8, (b) 1.0 (sc lattice) and (c) c/a = 1.2. The 
AFM propagation vector in the ordered AFM state at tA < 1 
is k = (i, |, 0) r.l.u. and the easy axis is the c axis [001] for 
both c/a = 0.8 and 1.2. The data were plotted using Eqs. (80) 
(Curie-Weiss law) for tA > 1, and (85) (xu) an d (87) (x||) for 
t A < 1. 



FIG. 15: (Color online) (a) Magnetic susceptibility x ver_ 
sus temperature T of tetragonal MnF 2 crystals for applied 
fields along the c axis ( Xc ) and in the ab plane ( Xab )- 62,63 
(b) Anisotropy Xc ~Xab versus T (solid blue squares). 64 Note 
the factor of 100 difference between the two ordinate scales 
in (a) and (b). The red solid curve is the MFT prediction for 
magnetic anisotropy arising from magnetic dipole interactions 
obtained using Eq. (96). 


to MFT, x_l = Xab for T < Xn should be independent of 
T, which is well satisfied. On the other hand, x\\ = Xc 
should go to zero as T — > 0, as also observed. We ob¬ 
tained a fairly good fit to X||(T < Tn) using our MFT 
with no adjustable parameters. 40 The fit function used 
was similar to the equation we obtained for X|| (T) for 
the pure dipole AFM in Eqs. (87) and Fig. 14. 

The anisotropy Ax(T) = Xc{T)—Xab{T) was measured 
with a torque magnetometer and the results are shown 
in Fig. 15(b). 64 The A\ data measured with the torque 
magnetometer 64 for T < Tn agree with the anisotropy 
calculated from the direct measurements 62,63 in Fig. 15. 
For T > Tn, a comparison of the data in Figs. 15(a) 
and 15(b) shows that |Ax|/x ~ 0.1% for T > Tn- From 
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Eq. (96), the anisotropy of x is predicted to be 

Axm(T) = ^^3 (Ao[ooi] — Ao[ioo])- (99) 

Using the values of a and Ao[ooil — ^o[ioo] i R Eqs. (97) 
and (98), respectively, and the xm(T) data in Fig. 15(a), 
A%m(T) was calculated from Eq. (99) and the result 
is shown as the solid red curve in Fig. 15(b) (see also 
Ref. 60). The calculation is in excellent agreement with 
the data for T > 150 K, suggesting that the MDI is re¬ 
sponsible for the x anisotropy in this T range, or at least 
reinforces this anisotropy. However, the data are increas¬ 
ingly suppressed to lower values below 130 K, which likely 
result from the onset of dynamic short-range collinear 
AFM correlations along the c axis with a correlation 
length that eventually diverges at Tn = 69 K, where 
from Fig. 15(a), A%m grows to become large and even 
more negative below that temperature. 

XII. SUMMARY 

A detailed summary of the paper is given in the Ab¬ 
stract to the paper. Here we provide a few additional 
comments. 

The eigenvalues and eigenvectors of the MDI tensor 
were determined for specified magnetic wavevectors and 
spin lattices. The eigenvalues give the energy of a spin in 
the magnetic fields of the local moments inside a Lorentz 
sphere of radius R in units of the a-axis lattice param¬ 
eter a. For 3D lattices, R/a = 50 was usually used, 
for a 2D circle R/a < 1000 and for a spin chain with 
R/a = oo the eigenvalues were determined exactly. The 
eigenvectors are the three orthogonal principal axis di¬ 
rections for collinear magnetic ordering. For uniaxial 3D 
spin lattices, these were calculated for c/a = 0.5 to 3 and 
the results presented in figures in the main text and Ap¬ 
pendix B and in tables in the Supplemental Material. 31 
We also calculated the eigenvalues and eigenvectors for 
noncollinear AFM structures including the 2D 120° tri¬ 
angular lattice and for the 2D and 3D coplanar non¬ 
collinear Shastry-Sutherland lattice and GdB 4 magnetic 
structure. We compared the ordering-direction predic¬ 
tions with data for some Mn +2 , (S = 5/2), Gd +3 and 
Eu +2 , (S = 7/2) compounds and found good agreement. 
Disagreement occurred for the itinerant AFM BaFe 2 As 2 
and for the coplanar noncollinear AFM GdB 4 , which in¬ 
dicates that a stronger anisotropy source must be present 
in these compounds that defeats the preferences of the 
MDI. 

A significant contribution of this paper was to ap¬ 
ply our formulation of the Weiss MFT 11,40 to predict 
many properties of the ordered and PM states arising 
from MDIs. These include the magnetic ordering tem¬ 
perature T m , the ordered moment, the magnetic heat 
capacity, and for AFMs the perpendicular critical field, 
the anisotropic magnetic susceptibility versus tempera¬ 
ture for T < Tn, and the parameters of the Curie-Weiss 


law for the anisotropic susceptibility for both FMs and 
AFMs at T > T m . Within MFT, the contributions of 
different molecular field sources to these properties are 
additive. This means that the same theory can be used 
to treat purely magnetic dipole magnets or spin systems 
containing both exchange and dipole interactions. We 
recently used the theory to separate the magnetic dipole 
and exchange contributions to the properties of the bet 
compound EuCu 2 Sb 2 with c/a = 2.4 and Tn = 5.1 K, 
which then allowed estimates of the Eu-Eu exchange in¬ 
teractions to be made. 29 
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Appendix A: Direct and Reciprocal Lattices 

1. Orthogonal Bravais Lattices 

In a Bravais spin lattice each spin position is a point 
of inversion symmetry with respect to the other spins. 
For orthogonal lattices which include as special cases the 
linear chain, the simple square lattice, sc, bcc, fee, simple 
tetragonal and bet lattices, the unit cell origins are at 

r 7 * 6 - c 

— = n a a H— rifcb H—n c c, (Ala) 

a a a 

where n a , rib and n c are positive or negative integers 
or 0. For all spin lattices, we normalize all spin positions 
and interspin distances by the a-axis lattice parameter a. 
For body-centered spin lattices one also has atoms at the 
body centers 

l = («. + 1) a+ (»s + 1) b + - a (n c + I) c, (Alb) 

where c/a = 1 for the bcc lattice. The central magnetic 
moment /q is placed at = 0 and hence the sum over 
neighbors jlj at positions Yj = r.j l in Eq. (16c) excludes 
the set ( n a ,rib,n c ) = (0,0,0) in Eq. (Ala). With our 
formulation, Gj(k) does not explicitly contain the lattice 
parameters a or c, and for tetragonal Bravais lattices just 
the dimensionless c/a ratio appears as in Eqs. (Al). 

The reciprocal-lattice vectors in reciprocal-lattice units 
are 

k = m ia* + 7712 b* + 7773 c*, (A2) 
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where the nii satisfy 0 < to, < 1 and the reciprocal- 
lattice translation vectors are 


* 2?r . 

a = — a, 

a 




and a, b and c are the corresponding direct-lattice trans¬ 
lation vectors. We normalize k by 1/a, yielding 


ka = 2tt[ m\k + TO 2 b 


c/a 


TO 3 C 


)■ 


(A4) 


Using Eqs. (A4) and (Ala), for the unit cell origins one 
has 


A 



FIG. 16: (Color online) In-plane hexagonal lattice translation 
unit vectors a and b of the direct lattice and a* and b* of the 
reciprocal lattice, respectively. 


k ■ Vji = 2n(rnin a + m 2 n b + m^ric) (A5a) 

and for the body-center positions 

n b + ” l " m3 ( n ° + 

( 

where the c/a ratio has canceled out of both expressions. 

The sum in Eq. (16c) gives an “extinction condition” 
for the contribution to the sum in Eq. (16c) of the body- 
centered spins in the bcc lattice in Eq. (A5b), where the 
contribution is zero if k-r^ is an odd multiple of ir/2 rad. 
This extinction occurs, for example, for AFM wavevec- 
tors 



kr ji = 2 tt 


toi( n a +- )+to 2 


where here the b axis is at an angle of 120° with respect 
to the positive x axis as shown in Fig. 16 and the rq 
are again positive or negative integers or zero. In two 
dimensions one sets n c = 0. In Cartesian coordinates the 
translation unit vectors are 

a = i, b = --i+—j, c = k. (All) 

A magnetic ordering wavevector k is written in terms 
of the respective simple hexagonal reciprocal lattice vec¬ 
tors as 


k = mia* + TO 2 b* + TO3C*, (Al2a) 


k = 






(A6) 


where the nii are chosen to satisfy 0 < rrii < 1 and the 
reciprocal lattice translation vectors are given by 


For such cases, according to Eq. (15) which assumes a 
collincar magnetic structure, the spins at the body cen¬ 
ters of the unit cells have zero ordered moment and they 
make no contribution to the dipolar interaction tensor in 
Eq. (16c). The interaction tensor is then the same as for 
a simple tetragonal lattice of moments with the same c/a 
ratio and k value. 

For the fee lattice the lattice points are at the positions 
in Eq. (Ala) and at 


r 

a 

r 

a 

r 

a 


1 


1 


n a + - a+( rib-f- b + O, 


n “4-)a + 0+ (n c + 


n b 


(A7) 

(AS) 

(A9) 


2tt 


47r 


27r - 


a* = — i+-=j , b* = —=j, c* = —k 


n/3' 


I a* I = |b*| = 


Vo 


(A12b) 


47T 
x/3 a 


= 1 a, 6-axis r.l.u., (Al2c) 


|c*| = — = 1 c-axis r.l.u. 
c 


(A12d) 


In terms of a* and b*, the direct lattice unit vectors are 


a = — |_(2a* — b*) b=^(2b*-a*), c = c*. 

v / 3 v 

(A13) 

The expression for k ■ r :ji is the same as in Eq. (A5a). 


with corresponding changes to the expressions for k ■ 


2. Simple Hexagonal (Triangular) Bravais Lattice 

The normalized vectors r ji for the simple hexagonal 
lattice with a = b are given by 

r ji - , c 1 c 
— = n a a + rib b H— n c c, 
a 


a 


(A10) 
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Appendix B: Figures Showing the Approach to the 
Large-Radius Asymptotic Eigenvalues for Magnetic 
Ordering on 2D and 3D Spin Lattices 



1/(R/a) 



FIG. 17: (a) Eigenvalue A( 0 ,o,o)[o,o,i] for FM spin alignment 
along the c axis versus the inverse of the circle radius R around 
the central moment in units of the square lattice parameter a 
for the 2D simple square lattice. The a and 6-axis eigenvalues 
are each equal to — A(o,o,o)[o,o,i] /2. (b) Deviation of the data 
from the fit. The “noise” is due to the discrete nature of 
the lattice, not to numerical inaccuracy. The lines in (b) are 
guides to the eye. 
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FIG. 18: (a) Eigenvalue A(i/2,i/2,o)[o,o,i] for the Neel-type 
AFM moment alignment along the c axis versus the circle 
radius R around the central moment in units of the square 
lattice parameter a for the 2D simple square lattice. The el¬ 
and 6-axis eigenvalues are each equal to —A(i/ 2 ,i/ 2 ,o)[o,o,il /2- 
(b) Deviation of the data from the fit. The “noise” is due to 
the discrete nature of the lattice, not to numerical inaccuracy. 
The lines in (b) are guides to the eye. 
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FIG. 19: Dependences of the eigenvalue A( 0 ,o,o)[o,o,i] on the 
inverse radius (R/a)^ 1 of the Lorentz sphere for FM moment 
alignments [k = (0,0,0)] along the c axis in 3D simple tetrag¬ 
onal spin lattices with (a) c/a = 1.5 and (b) c/a = 3. The 
lines are guides to the eye. 



FIG. 20: Dependences of the eigenvalue A( 1 / 2 ,i/2,i/2)[o,o,i] for 
Neel-type ordering with k = (1/2,1/2,1/2) on the inverse ra¬ 
dius ( R/a ) -1 of the Lorentz sphere for AFM moment align¬ 
ments along the c axis in 3D simple tetragonal spin lattices 
with (a) c/a = 1.5 and (b) c/a = 3. The lines are guides to 
the eye. The averages for R/a = 51-100 are shown. With in¬ 
creasing c/a, the averages of A(i/ 2 ,i/2,i/2)[o,o,i] for R/a = 51- 
100 approach the 2D square-lattice limit A(i/ 2 ] i/ 2 ,o)[o,o,i] = 
2.645 887 in Table I, as shown in Fig. 3(b). 
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Appendix C: Figures Showing Dipolar Eigenvectors 
and Eigenvalues versus the c/a Ratio for Tetragonal 
and Hexagonal Bravais Spin Lattices and for the 
Honeycomb Lattice 



c/a 

FIG. 21: (Color online) Dependence of the eigenvalue 
A(o,o,o)[o,o,i] on the c/a ratio for a simple tetragonal lattice 
with a FM alignment of the magnetic moments along the 
c axis. From the figure, one sees that FM alignment along 
the c axis is the most stable for c/a < 1, but for c/a > 1 FM 
alignment along the a or b axis is energetically favorable. 




FIG. 23: (Color online) Eigenvalues (a) Ap/ 2 , 1 / 2 , 0 ) f° r AFM 
wavevector k = (1/2,1/2,0) r.l.u. and (b) A( 0 j o,i/ 2 ) f° r AFM 
wavevector k = (0,0,1/2) r.l.u. versus the c/a ratio for a simple 
tetragonal lattice with the moments aligned along [1, 0, 0] or 
[0, 1, 0] (a or b axis, solid blue squares) or [0, 0, 1] (c axis, 
solid green diamonds). 


FIG. 22: (Color online) Eigenvalues for wavevector k = 
(1/2,0,0) r.l.u. versus the c/a ratio for a simple tetragonal 
or bet lattice with the moments aligned along [010] (6 axis, 
solid red circles), [001] (c axis, solid green diamonds) or [100] 
(a axis, solid blue squares). The 2D limits for c/a —> 00 are 
shown as horizontal dashed lines. 
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c/a 

FIG. 24: (Color online) Eigenvalues for wavevector k = 
(1/2,1/2,0) r.l.u. versus the c/a ratio for a bet spin lattice with 
the moments aligned along [1, —1, 0] (solid red circles), [001] 
(c axis, solid green diamonds) or [110] (solid blue squares). 



FIG. 25: (Color online) Eigenvalues for wavevectors (a) k = 
0 (FM) or (1,0,0) and (b) k = (1/2,1/2,0) r.l.u. versus the c/a 
ratio for a simple hexagonal (stacked triangular) spin lattice 
with the moments aligned along the indicated principal axes. 
The 2D limits of the respective eigenvalues for c/a —» oo are 
shown by horizontal dashed lines. 
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FIG. 26: (Color online) Eigenvalues for wavevectors (a) k = 
(1/3,1/3,1/3) and (b) k = (1/2,1/2,1/2) r.l.u. versus the c/a 
ratio for a simple hexagonal (stacked triangular) spin lattice 
with the moments aligned along the indicated principal axes. 
The 2D limits of the respective eigenvalues for c/a —» oo are 
shown by horizontal dashed lines. 



FIG. 27: (Color online) Eigenvalues for wavevectors (a) k = 
(1/3,1/3,0) and (b) k = (1/3,1/3,1/2) r.l.u. versus the c/a 
ratio for a simple hexagonal (stacked triangular) spin lattice 
with the moments aligned along the indicated principal axes. 
The 2D limits of the respective eigenvalues for c/a —» oo are 
shown by horizontal dashed lines. 
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c/a 



FIG. 28: (Color online) Eigenvalues for propagation vectors 
(a) k = (0,0,0) (FM) and (b) k = (1/2,0,0) r.l.u. versus the c/a 
ratio for a honeycomb spin lattice with the moments aligned 
along the indicated principal axes. The 2D limits of the re¬ 
spective eigenvalues for c/a —> oo are shown by horizontal 
dashed lines. 



FIG. 29: (Color online) Eigenvalues for AFM propagation 
vectors (a) k = (0,0,1/2) (Neel-type in all directions) and (b) 
k = (0,0,0) r.l.u. (Neel-type in ab plane and FM-alignment 
along c axis) versus the c/a ratio for a honeycomb spin lattice 
with the moments aligned along the indicated principal axes. 
The 2D limits of the respective eigenvalues for c/a —» oo are 
shown by horizontal dashed lines. 
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FIG. 30: (Color online) Eigenvalues for AFM propagation 
vectors k = (0,0,1/2) (FM alignment intraplane and AFM 
alignment interplane) versus the c/a ratio for a honeycomb 
spin lattice with the moments aligned along the indicated 
principal axes. The 2D limits of the respective eigenvalues 
for c/a —> oo are shown by horizontal dashed lines. 
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Supplementary Information. Tables of Dipolar 
Eigenvectors and Eigenvalues versus the c/a Ratio 
for Tetragaonal and Hexagonal Bravais Spin Lattices 
and for the Non-Bravais Honeycomb Lattice 


TABLE VII: Simple Tetragonal Spin Lattices with Fer¬ 
romagnetic Alignment and k = (0,0,0). Eigenvalues Aka 
and eigenvectors jl = [fj, x , fi y , y z ] in Cartesian coordinates 
of the magnetic dipole interaction tensor Gi(k) in Eq. (16c) 
for various values of the magnetic wavevector k in reciprocal 
lattice units (r.l.u.) for simple tetragonal spin lattices with 
collinear magnetic moment alignments. The most positive 
Aka value(s) corresponds to the lowest energy value accord¬ 
ing to Eq. (16d). The Cartesian x, y and z axes are along 
the a, b and c axes of the tetragonal lattices, respectively. 
The accuracy of the values is estimated to be < ±0.001. Also 
shown are the differences between the eigenvalues for different 
ordering axes for a given k, which determine the anisotropy 
energies via Eq. (16d). 


c/a 

[001] 

\o,0, 0 )a 
[100], [010] 

[100] - [001] 

0.5 

30.0834 

-15.0417 

-45.1250 

0.6 

15.2488 

-7.6244 

-22.8732 

0.7 

7.9376 

-3.9688 

-11.9064 

0.8 

3.9382 

-1.9691 

-5.9073 

0.9 

1.5482 

-0.7741 

-2.3222 

1 

0 

0 

0 

1.1 

-1.0757 

0.5378 

1.6135 

1.2 

-1.8728 

0.9364 

2.8091 

1.3 

-2.4947 

1.2474 

3.7421 

1.4 

-2.9994 

1.4997 

4.4991 

1.5 

-3.4217 

1.7109 

5.1326 

1.6 

-3.7830 

1.8915 

5.6745 

1.7 

-4.0975 

2.0487 

6.1462 

1.8 

-4.3750 

2.1875 

6.5625 

1.9 

-4.6214 

2.3107 

6.9321 

2 

-4.8436 

2.4218 

7.2653 

2.1 

-5.0434 

2.5217 

7.5651 

2.2 

-5.2249 

2.6125 

7.8374 

2.3 

-5.3900 

2.6950 

8.0850 

2.4 

-5.5425 

2.7712 

8.3137 

2.5 

-5.6823 

2.8411 

8.5234 

2.6 

-5.8110 

2.9055 

8.7165 

2.7 

-5.9303 

2.9652 

8.8955 

2.8 

-6.0413 

3.0207 

9.0620 

2.9 

-6.1447 

3.0724 

9.2171 

3 

-6.2410 

3.1205 

9.3616 
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TABLE VIII: Simple Tetragonal Antiferromagnetic Spin Lattices with k = (1/2,0,1/2) and (1/2,1/2,0) r.l.u. The 

symbol descriptions are the same as in the caption to Table VII. 


Aka 

c/a 

(l,o,|) 

[100] 

[010] 

[001] 

[010]-[001] 

[010]-[100] 

(U,o) 
[100], [010] 

[001] 

[001]-[100] 

0.5 

13.6119 

15.1664 

-28.7783 

43.9447 

1.5545 

-19.236 

38.4719 

57.7079 

0.6 

6.7759 

9.7270 

-16.5029 

26.2299 

2.9511 

-11.1459 

22.2918 

33.4377 

0.7 

2.8478 

7.2967 

-10.1445 

17.4412 

4.4489 

-7.0563 

14.1126 

21.1689 

0.8 

0.3134 

6.1612 

-6.4746 

12.6358 

5.8478 

-4.7977 

9.5955 

14.3932 

0.9 

-1.4328 

5.6177 

-4.1849 

9.8026 

7.0505 

-3.4771 

6.9543 

10.4314 

1 

-2.6767 

5.3535 

-2.6768 

8.0303 

8.0302 

-2.6767 

5.3535 

8.0302 

1.1 

-3.5784 

5.2246 

-1.6462 

6.8708 

8.8030 

-2.1807 

4.3614 

6.5421 

1.2 

-4.2360 

5.1611 

-0.9251 

6.0862 

9.3971 

-1.8691 

3.7381 

5.6072 

1.3 

-4.7177 

5.1295 

-0.4118 

5.5413 

9.8472 

-1.6716 

3.3431 

5.0147 

1.4 

-5.0709 

5.1143 

-0.0434 

5.1577 

10.1852 

-1.5462 

3.0925 

4.6388 

1.5 

-5.3294 

5.1064 

0.2230 

4.8834 

10.4358 

-1.4655 

2.9311 

4.3966 

1.6 

-5.5185 

5.1026 

0.4159 

4.6867 

10.6211 

-1.4143 

2.8286 

4.2429 

1.7 

-5.6572 

5.1002 

0.5570 

4.5432 

10.7574 

-1.3815 

2.7631 

4.1446 

1.8 

-5.7585 

5.1001 

0.6584 

4.4417 

10.8586 

-1.3605 

2.7209 

4.0814 

1.9 

-5.8328 

5.0996 

0.7332 

4.3664 

10.9324 

-1.3473 

2.6946 

4.0419 

2 

-5.8871 

5.0991 

0.7880 

4.3111 

10.9862 

-1.3384 

2.6767 

4.0150 

2.1 

-5.9269 

5.0993 

0.8276 

4.2717 

11.0262 

-1.3328 

2.6656 

3.9984 

2.2 

-5.9559 

5.0992 

0.8567 

4.2425 

11.0551 

-1.3292 

2.6585 

3.9878 

2.3 

-5.9771 

5.0991 

0.8780 

4.2211 

11.0762 

-1.3269 

2.6537 

3.9806 

2.4 

-5.9924 

5.0992 

0.8932 

4.2060 

11.0916 

-1.3255 

2.6511 

3.9766 

2.5 

-6.0038 

5.0992 

0.9046 

4.1946 

11.1030 

-1.3247 

2.6494 

3.9741 

2.6 

-6.0118 

5.0988 

0.9130 

4.1858 

11.1106 

-1.3239 

2.6478 

3.9717 

2.7 

-6.0180 

5.0989 

0.9191 

4.1798 

11.1169 

-1.3235 

2.6471 

3.9706 

2.8 

-6.0222 

5.0990 

0.9232 

4.1758 

11.1212 

-1.3235 

2.6471 

3.9706 

2.9 

-6.0256 

5.0988 

0.9268 

4.1720 

11.1244 

-1.3234 

2.6469 

3.9703 

3 

-6.0278 

5.0989 

0.9289 

4.1700 

11.1267 

-1.3231 

2.6461 

3.9692 
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TABLE IX: Simple Tetragonal and Body-Centered 
Tetragonal Spin Lattice with k = (0,0,1/2) r.l.u. The 

symbol descriptions are the same as in the caption to Ta¬ 
ble VII. 


c/a 

[100], [010] 

[001] 

[100]-[001] 

0.5 

15.0417 

-30.0834 

45.1251 

0.6 

9.4972 

-18.9943 

28.4915 

0.7 

6.9682 

-13.9364 

20.9046 

0.8 

5.7521 

-11.5042 

17.2563 

0.9 

5.1491 

-10.2982 

15.4473 

1 

4.8436 

-9.6871 

14.5307 

1.1 

4.6871 

-9.3742 

14.0613 

1.2 

4.6061 

-9.2122 

13.8183 

1.3 

4.5637 

-9.1273 

13.6910 

1.4 

4.5416 

-9.0833 

13.6249 

1.5 

4.5302 

-9.0604 

13.5906 

1.6 

4.5238 

-9.0476 

13.5714 

1.7 

4.5205 

-9.0411 

13.5616 

1.8 

4.5187 

-9.0373 

13.5560 

1.9 

4.5184 

-9.0367 

13.5550 

2 

4.5173 

-9.0346 

13.5519 

2.1 

4.5174 

-9.0347 

13.5521 

2.2 

4.5168 

-9.0337 

13.5506 

2.3 

4.5171 

-9.0343 

13.5514 

2.4 

4.5163 

-9.0325 

13.5487 

2.5 

4.5167 

-9.0333 

13.5500 

2.6 

4.5168 

-9.0337 

13.5506 

2.7 

4.5169 

-9.0338 

13.5507 

2.8 

4.5167 

-9.0333 

13.5500 

2.9 

4.5168 

-9.0336 

13.5504 

3 

4.5175 

-9.0351 

13.5527 
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TABLE X: Simple Tetragonal and Body-Centered Tetragonal Spin Lattices with k = (1/2,0,0) and 
(1/2,1/2,1/2) r.l.u. The symbol descriptions are the same as in the caption to Table VII. 


•^kct 

c/a 

(3.0,0) 

[100] 

[010] 

[001] 

[010]-[001] 

[010] —[100] 

(I I I) 

V 2 ’ 2’ 2) 
[100], [010] 

[001] 

[001]-[100] 

0.5 

-32.9866 

-5.4794 

38.4660 

-43.9454 

27.5072 

13.8825 

-27.7650 

-41.6475 

0.6 

-22.6079 

0.3470 

22.2609 

-21.9139 

22.9549 

7.4078 

-14.8156 

-22.2234 

0.7 

-16.8911 

2.8706 

14.0205 

-11.1499 

19.7617 

3.9638 

-7.9276 

-11.8914 

0.8 

-13.4266 

4.0289 

9.3977 

-5.3688 

17.4555 

1.9691 

-3.9383 

-5.9074 

0.9 

-11.1910 

4.5782 

6.6128 

-2.0346 

15.7692 

0.7562 

-1.5124 

-2.2686 

1 

-9.6874 

4.8437 

4.8437 

0.0000 

14.5311 

0.0000 

0.0000 

0.0000 

1.1 

-8.6477 

4.9729 

3.6748 

1.2981 

13.6206 

-0.4778 

0.9556 

1.4334 

1.2 

-7.9160 

5.0370 

2.8790 

2.1580 

12.9530 

-0.7821 

1.5641 

2.3461 

1.3 

-7.3954 

5.0683 

2.3271 

2.7412 

12.4637 

-0.9763 

1.9526 

2.9289 

1.4 

-7.0217 

5.0837 

1.9380 

3.1457 

12.1054 

-1.1009 

2.2018 

3.3027 

1.5 

-6.7521 

5.0915 

1.6606 

3.4309 

11.8436 

-1.1805 

2.3611 

3.5416 

1.6 

-6.5566 

5.0950 

1.4616 

3.6334 

11.6516 

-1.2315 

2.4630 

3.6945 

1.7 

-6.4150 

5.0969 

1.3181 

3.7788 

11.5119 

-1.2644 

2.5287 

3.7931 

1.8 

-6.3119 

5.0980 

1.2139 

3.8841 

11.4099 

-1.2852 

2.5704 

3.8556 

1.9 

-6.2367 

5.0986 

1.1381 

3.9605 

11.3353 

-1.2987 

2.5974 

3.8961 

2 

-6.1821 

5.0986 

1.0835 

4.0151 

11.2807 

-1.3074 

2.6147 

3.9220 

2.1 

-6.1421 

5.0985 

1.0436 

4.0549 

11.2406 

-1.3132 

2.6263 

3.9395 

2.2 

-6.1131 

5.0988 

1.0143 

4.0845 

11.2119 

-1.3169 

2.6339 

3.9509 

2.3 

-6.0917 

5.0988 

0.9929 

4.1059 

11.1905 

-1.3188 

2.6376 

3.9564 

2.4 

-6.0760 

5.0990 

0.9770 

4.1220 

11.1750 

-1.3201 

2.6402 

3.9603 

2.5 

-6.0652 

5.0992 

0.9660 

4.1332 

11.1644 

-1.3213 

2.6426 

3.9639 

2.6 

-6.0566 

5.0989 

0.9577 

4.1412 

11.1555 

-1.3218 

2.6437 

3.9655 

2.7 

-6.0507 

5.0988 

0.9519 

4.1469 

11.1495 

-1.3222 

2.6444 

3.9666 

2.8 

-6.0463 

5.0990 

0.9473 

4.1517 

11.1453 

-1.3223 

2.6445 

3.9668 

2.9 

-6.0430 

5.0988 

0.9442 

4.1546 

11.1418 

-1.3226 

2.6452 

3.9678 

3 

-6.0408 

5.0990 

0.9418 

4.1572 

11.1398 

-1.3227 

2.6455 

3.9682 
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TABLE XI: Body-Centered Tetragonal Spin Lattice 
with ferromagnetic moment alignment and k = 
(0,0,0) r.l.u. The symbol descriptions are the same as in 
the caption to Table VII. 


c/a 

[001] 

-^ka 

[100], [010] 

[100]-[001] 

0.50 

21.9948 

-10.9974 

-32.9923 

0.60 

9.0718 

-4.5359 

-13.6077 

0.70 

3.5213 

-1.7607 

-5.2820 

0.80 

1.1666 

-0.5833 

-1.7499 

0.90 

0.2635 

-0.1317 

-0.3952 

0.93 

0.1405 

-0.0702 

-0.2107 

0.95 

0.0825 

-0.0412 

-0.1237 

0.97 

0.0404 

-0.0202 

-0.0606 

0.99 

0.0105 

-0.0052 

-0.0157 

1 

0 

0 

0 

1.01 

-0.0083 

0.0042 

0.0125 

1.03 

-0.0182 

0.0091 

0.0273 

1.05 

-0.0231 

0.0115 

0.0346 

1.07 

-0.0219 

0.0109 

0.0328 

1.09 

-0.0170 

0.0085 

0.0255 

1.10 

-0.0143 

0.0071 

0.0214 

1.11 

-0.0098 

0.0049 

0.0147 

1.13 

-0.0015 

0.0007 

0.0022 

1.15 

0.0084 

-0.0042 

-0.0127 

1.17 

0.0179 

-0.0089 

-0.0268 

1.19 

0.0269 

-0.0135 

-0.0404 

1.20 

0.0307 

-0.0154 

-0.0461 

1.21 

0.0346 

-0.0173 

-0.0520 

1.23 

0.0413 

-0.0206 

-0.0619 

1.25 

0.0463 

-0.0232 

-0.0695 

1.27 

0.0493 

-0.0246 

-0.0739 

1.29 

0.0497 

-0.0248 

-0.0745 

1.30 

0.0491 

-0.0246 

-0.0737 

1.31 

0.0489 

-0.0244 

-0.0733 

1.33 

0.0448 

-0.0224 

-0.0672 

1.35 

0.0379 

-0.0190 

-0.0569 

1.37 

0.0294 

-0.0147 

-0.0440 

1.39 

0.0175 

-0.0088 

-0.0263 

1.40 

0.0110 

-0.0055 

-0.0165 

1.41 

0.0039 

-0.0019 

-0.0058 

V2 

0 

0 

0 

1.43 

-0.0122 

0.0061 

0.0183 

1.45 

-0.0324 

0.0162 

0.0486 

1.47 

-0.0535 

0.0268 

0.0803 

1.49 

-0.0764 

0.0382 

0.1146 

1.50 

-0.0896 

0.0448 

0.1344 

1.60 

-0.2452 

0.1226 

0.3678 

1.70 

-0.4435 

0.2217 

0.6652 

1.80 

-0.6731 

0.3366 

1.0097 

1.90 

-0.9220 

0.4610 

1.3830 

2.00 

-1.1823 

0.5911 

1.7734 

2.10 

-1.4452 

0.7226 

2.1678 

2.20 

-1.7061 

0.8530 

2.5591 

2.30 

-1.9611 

0.9806 

2.9417 

2.40 

-2.2097 

1.1048 

3.3145 

2.50 

-2.4475 

1.2238 

3.6713 

2.60 

-2.6739 

1.3369 

4.0108 

2.70 

-2.8900 

1.4450 

4.3349 

2.80 

-3.0954 

1.5477 

4.6431 

2.90 

-3.2898 

1.6449 

4.9348 

3.00 

-3.4735 

1.7367 

5.2102 
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TABLE XII: Body-Centered Tetragonal Spin Lattices with k = (1/2,0,1/2) r.l.u. Eigenvalues Aka and eigenvectors 


jj, = [p X! fi y , fi z ] in Cartesian coordinates of the magnetic dipole interaction tensor Gi(k) in Eq. (16c) with collinear magnetic 
moment alignments. The accuracy of the eigenvalues is estimated to be < ±0.001. Also shown are two differences between the 
eigenvalues for the three different ordering axes. 


A 

c/a 

X 

[010] 

[y/1 — X 2 ,0, —x] 

[: X , 0, \/l — X 2 ] 

[010] — [VI — * 2 , 0, -x] 

[ 010 ] - {x,0,Vl~x 2 } 

0.5 

0.15297 

15.1665 

14.6525 

-29.8189 

0.5140 

44.9854 

0.6 

0.32332 

9.7270 

9.8527 

-19.5797 

-0.1257 

29.3067 

0.7 

0.48513 

7.2967 

8.6251 

-15.9218 

-1.3284 

23.2185 

0.8 

0.59407 

6.1612 

8.4577 

-14.6189 

-2.2965 

20.7801 

0.9 

0.66190 

5.6177 

8.3111 

-13.9288 

-2.6934 

19.5465 

1 

0.70711 

5.3535 

7.9437 

-13.2971 

-2.5902 

18.6506 

1.1 

0.74051 

5.2246 

7.3770 

-12.6016 

-2.1524 

17.8262 

1.2 

0.76766 

5.1611 

6.6884 

-11.8495 

-1.5273 

17.0106 

1.3 

0.79145 

5.1295 

5.9492 

-11.0787 

-0.8197 

16.2082 

1.4 

0.81346 

5.1143 

5.2144 

-10.3287 

-0.1001 

15.4430 

V2 

0.81650 

5.1127 

5.1127 

-10.2254 

0 

15.3381 

1.5 

0.83442 

5.1064 

4.5196 

-9.6260 

0.5868 

14.7324 

1.6 

0.85465 

5.1026 

3.8873 

-8.9899 

1.2153 

14.0925 

1.7 

0.87419 

5.1002 

3.3301 

-8.4303 

1.7701 

13.5305 

1.8 

0.89285 

5.1001 

2.8480 

-7.9481 

2.2521 

13.0482 

1.9 

0.91035 

5.0996 

2.4435 

-7.5431 

2.6561 

12.6427 

2 

0.92634 

5.0991 

2.1105 

-7.2096 

2.9886 

12.3087 

2.1 

0.94057 

5.0993 

1.8402 

-6.9395 

3.2591 

12.0388 

2.2 

0.95284 

5.0992 

1.6258 

-6.7250 

3.4734 

11.8242 

2.3 

0.96316 

5.0991 

1.4576 

-6.5567 

3.6415 

11.6558 

2.4 

0.97160 

5.0992 

1.3273 

-6.4266 

3.7719 

11.5258 

2.5 

0.97836 

5.0992 

1.2281 

-6.3273 

3.8711 

11.4265 

2.6 

0.98367 

5.0988 

1.1528 

-6.2517 

3.9460 

11.3505 

2.7 

0.98776 

5.0989 

1.0965 

-6.1954 

4.0024 

11.2943 

2.8 

0.99088 

5.0990 

1.0540 

-6.1530 

4.0450 

11.2520 

2.9 

0.99324 

5.0988 

1.0231 

-6.1220 

4.0757 

11.2208 

3 

0.99501 

5.0989 

0.9996 

-6.0985 

4.0993 

11.1974 
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TABLE XIII: Body-Centered Tetragonal Spin Lattices with k= (1/2,1/2,0) and (0,0,1) r.l.u. Eigenvalues Aka and 
eigenvectors fi = [n x , fi y , )i z ] in Cartesian coordinates of the magnetic dipole interaction tensor Gj(k) in Eq. (16c) for two values 
of the magnetic wavevector k in reciprocal lattice units (r.l.u.) and collinear magnetic moment alignments. The accuracy of 
the eigenvalues is estimated to be < ±0.001. Also shown are the differences between the eigenvalues for different ordering axes 
for a given k. 


c/a 

(il.O) 

[001] 

[no] 

[110] 

[001]—[110] 

[001]—[110] 

(0,0,1) 
[100],[010] 

[001] 

[100]-[001] 

0.5 

38.4705 

8.0868 

-46.5573 

30.3837 

85.0278 

-19.086 

38.172 

-57.258 

0.6 

22.2918 

11.2533 

-33.5451 

11.0385 

55.8369 

-10.713 

21.426 

-32.139 

0.7 

14.1126 

11.5369 

-25.6495 

2.5757 

39.7621 

-6.177 

12.354 

-18.531 

0.8 

9.5955 

10.6756 

-20.2711 

-1.0801 

29.8666 

-3.355 

6.710 

-10.065 

0.9 

6.9543 

9.3705 

-16.3248 

-2.4162 

23.2791 

-1.416 

2.833 

-4.249 

1 

5.3534 

7.9437 

-13.2971 

-2.5903 

18.6505 

0.000 

0.000 

0.000 

1.1 

4.3614 

6.5552 

-10.9166 

-2.1938 

15.2780 

1.069 

-2.137 

3.206 

1.2 

3.7381 

5.2829 

-9.0210 

-1.5448 

12.7591 

1.888 

-3.776 

5.664 

1.3 

3.3431 

4.1586 

-7.5017 

-0.8155 

10.8448 

2.519 

-5.039 

7.558 

1.4 

3.0925 

3.1897 

-6.2822 

-0.0972 

9.3747 

3.005 

-6.010 

9.015 

V2 

3.0640 

3.0640 

-6.1280 

0 

9.1920 

3.064 

-6.128 

9.192 

1.5 

2.9311 

2.3693 

-5.3004 

0.5618 

8.2315 

3.377 

-6.754 

10.131 

1.6 

2.8286 

1.6843 

-4.5129 

1.1443 

7.3415 

3.660 

-7.321 

10.981 

1.7 

2.7627 

1.1176 

-3.8803 

1.6451 

6.6430 

3.876 

-7.752 

11.627 

1.8 

2.7209 

0.6513 

-3.3722 

2.0696 

6.0931 

4.038 

-8.077 

12.115 

1.9 

2.6946 

0.2700 

-2.9646 

2.4246 

5.6592 

4.160 

-8.321 

12.481 

2 

2.6765 

-0.0384 

-2.6381 

2.7149 

5.3146 

4.252 

-8.505 

12.757 

2.1 

2.6656 

-0.2896 

-2.3760 

2.9552 

5.0416 

4.321 

-8.642 

12.962 

2.2 

2.6585 

-0.4921 

-2.1664 

3.1506 

4.8249 

4.372 

-8.744 

13.116 

2.3 

2.6537 

-0.6555 

-1.9982 

3.3092 

4.6519 

4.409 

-8.819 

13.228 

2.4 

2.6511 

-0.7876 

-1.8635 

3.4387 

4.5146 

4.438 

-8.875 

13.313 

2.5 

2.6494 

-0.8938 

-1.7556 

3.5432 

4.4050 

4.459 

-8.917 

13.376 

2.6 

2.6478 

-0.9783 

-1.6695 

3.6261 

4.3173 

4.474 

-8.948 

13.422 

2.7 

2.6471 

-1.0462 

-1.6009 

3.6933 

4.2480 

4.485 

-8.971 

13.456 

2.8 

2.6471 

-1.1014 

-1.5457 

3.7485 

4.1928 

4.494 

-8.987 

13.481 

2.9 

2.6469 

-1.1461 

-1.5008 

3.7930 

4.1477 

4.500 

-9.000 

13.499 

3 

2.6461 

-1.1803 

-1.4658 

3.8264 

4.1119 

4.504 

-9.009 

13.513 
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TABLE XIV: Simple 3D Hexagonal (Stacked Triangular) Spin Lattices with k = (1/2,1/2,0) and (0,0,0) or 
(1,0,0) r.l.u. Eigenvalues Aka and eigenvectors /t = Hy, Hz\ in Cartesian coordinates of the magnetic dipole interaction 
tensor G, (k) in Eq. (16c) for two values of the magnetic wavevector k in reciprocal lattice units (r.l.u.) and collinear magnetic 
moment alignments. The x axis is in the direction of the hexagonal a axis. The accuracy of the eigenvalues is estimated to be 
< ±0.001. Also shown are the differences between the eigenvalues for different ordering axes for a given k. 


Aka 

c/a 

(U.O) 

[ 001 ] 

“12a 

1 

O 

[ 3 .^. 0 ] 

[ 3 .^. 0 ] 

-[001] 

1 ±3 0 ] 

2 ’ 2 > U J 

'±1 _1 q| 

2 1 2> u | 

( 0 , 0 , 0 ), ( 1 , 0 , 0 ) 

[ 001 ] 

[ 100 ], [ 010 ] 

[ 100 ] - [ 001 ] 

0.5 

38.4698 

-31.0341 

-7.4357 

-45.9055 

23.5983 

28.7829 

-14.3914 

-43.1743 

0.6 

22.2769 

-20.9871 

-1.2899 

-23.5668 

19.6972 

14.1513 

-7.0757 

-21.2270 

0.7 

14.0669 

-15.5182 

1.4513 

-12.6156 

16.9695 

6.9640 

-3.4820 

-10.4460 

0.8 

9.4972 

-12.2545 

2.7573 

-6.7399 

15.0118 

3.0273 

-1.5137 

-4.5410 

0.9 

6.7848 

-10.1901 

3.4053 

-3.3795 

13.5953 

0.6542 

-0.3271 

-0.9813 

1 

5.1009 

-8.8361 

3.7353 

-1.3656 

12.5714 

-0.9104 

0.4552 

1.3656 

1.1 

4.0221 

-7.9278 

3.9056 

-0.1165 

11.8334 

-2.0227 

1.0113 

3.0340 

1.2 

3.3153 

-7.3096 

3.9943 

0.6789 

11.3038 

-2.8679 

1.4339 

4.3018 

1.3 

2.8443 

-6.8859 

4.0416 

1.1972 

10.9275 

-3.5418 

1.7709 

5.3128 

1.4 

2.5278 

-6.5942 

4.0664 

1.5386 

10.6606 

-4.1002 

2.0501 

6.1503 

1.5 

2.3125 

-6.3921 

4.0797 

1.7672 

10.4718 

-4.5732 

2.2866 

6.8599 

1.6 

2.1656 

-6.2521 

4.0865 

1.9209 

10.3386 

-4.9823 

2.4911 

7.4734 

1.7 

2.0650 

-6.1551 

4.0901 

2.0251 

10.2452 

-5.3406 

2.6703 

8.0110 

1.8 

1.9952 

-6.0874 

4.0922 

2.0970 

10.1796 

-5.6583 

2.8291 

8.4874 

1.9 

1.9476 

-6.0410 

4.0934 

2.1458 

10.1343 

-5.9423 

2.9712 

8.9135 

2 

1.9144 

-6.0082 

4.0938 

2.1794 

10.1019 

-6.1973 

3.0986 

9.2959 

2.1 

1.8915 

-5.9856 

4.0941 

2.2026 

10.0798 

-6.4280 

3.2140 

9.6420 

2.2 

1.8757 

-5.9700 

4.0944 

2.2187 

10.0644 

-6.6369 

3.3184 

9.9553 

2.3 

1.8641 

-5.9587 

4.0946 

2.2305 

10.0533 

-6.8282 

3.4141 

10.2423 

2.4 

1.8569 

-5.9512 

4.0943 

2.2375 

10.0456 

-7.0036 

3.5018 

10.5055 

2.5 

1.8515 

-5.9458 

4.0943 

2.2428 

10.0402 

-7.1647 

3.5824 

10.7471 

2.6 

1.8477 

-5.9421 

4.0944 

2.2467 

10.0365 

-7.3135 

3.6567 

10.9702 

2.7 

1.8452 

-5.9396 

4.0944 

2.2493 

10.0340 

-7.4513 

3.7256 

11.1769 

2.8 

1.8432 

-5.9378 

4.0946 

2.2513 

10.0324 

-7.5797 

3.7898 

11.3695 

2.9 

1.8419 

-5.9366 

4.0947 

2.2528 

10.0313 

-7.6983 

3.8491 

11.5474 

3 

1.8410 

-5.9360 

4.0950 

2.2540 

10.0311 

-7.8098 

3.9049 

11.7147 
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TABLE XV: Simple 3D Hexagonal (Stacked Triangular) Spin Lattices with k = (1/2,1/2,1/2) and (1/3,1/3,1/3). 

Eigenvalues Aka and eigenvectors fi = [fi x , fi y , fi z \ in Cartesian coordinates of the magnetic dipole interaction tensor Gi(k) in 
Eq. (16c) for two values of the magnetic wavevector k in reciprocal lattice units (r.l.u.) and collinear magnetic moment 
alignments. The x axis is in the direction of the hexagonal a axis. The accuracy of the eigenvalues is estimated to be < ±0.001. 
Also shown are the differences between the eigenvalues for different ordering axes for a given k. 


•^kcr 

c/a 

(i I I) 

1 2 ’ 2’ 2/ 

[001] 



-[001] 

i ^ ol 

2’ 2 ’ U J 

±1 _I ol 

2 ’ 2’ v \ 

(1 I I) 

1 3 ’ 3 5 3/ 

[001] 

[100], [010] 

[001] - [100] 

0.5 

-28.2681 

14.9063 

13.3618 

43.1743 

1.5444 

-15.2234 

7.6117 

-22.8352 

0.6 

-15.6463 

9.2784 

6.3679 

24.9247 

2.9105 

-7.4986 

3.7493 

-11.2479 

0.7 

-9.0113 

6.6815 

2.3298 

15.6928 

4.3517 

-3.5057 

1.7528 

-5.2585 

0.8 

-5.1660 

5.4157 

-0.2497 

10.5817 

5.6653 

-1.2648 

0.6324 

-1.8972 

0.9 

-2.7943 

4.7791 

-1.9848 

7.57341 

6.7638 

0.0638 

-0.0319 

0.0957 

1 

-1.2721 

4.4523 

-3.1802 

5.72436 

7.6325 

0.8795 

-0.4397 

1.3192 

1.1 

-0.2708 

4.2831 

-4.0123 

4.55388 

8.2954 

1.3929 

-0.6965 

2.0894 

1.2 

0.3984 

4.1944 

-4.5928 

3.79602 

8.7873 

1.7210 

-0.8605 

2.5816 

1.3 

0.8511 

4.1472 

-4.9984 

3.29611 

9.1456 

1.9333 

-0.9666 

2.8999 

1.4 

1.1587 

4.1227 

-5.2814 

2.96403 

9.4042 

2.0710 

-1.0355 

3.1065 

1.5 

1.3696 

4.1103 

-5.4799 

2.74075 

9.5903 

2.1607 

-1.0803 

3.2410 

1.6 

1.5145 

4.1027 

-5.6172 

2.58825 

9.7199 

2.2194 

-1.1097 

3.3291 

1.7 

1.6147 

4.0985 

-5.7132 

2.48379 

9.8116 

2.2580 

-1.1290 

3.3870 

1.8 

1.6835 

4.0972 

-5.7807 

2.41363 

9.8779 

2.2834 

-1.1417 

3.4252 

1.9 

1.7305 

4.0971 

-5.8276 

2.36665 

9.9247 

2.2998 

-1.1499 

3.4496 

2 

1.7641 

4.0958 

-5.8599 

2.33175 

9.9557 

2.3108 

-1.1554 

3.4662 

2.1 

1.7870 

4.0953 

-5.8823 

2.30833 

9.9776 

2.3182 

-1.1591 

3.4773 

2.2 

1.8030 

4.0947 

-5.8977 

2.29180 

9.9924 

2.3230 

-1.1615 

3.4845 

2.3 

1.8141 

4.0945 

-5.9086 

2.28032 

10.0031 

2.3261 

-1.1631 

3.4892 

2.4 

1.8219 

4.0946 

-5.9165 

2.27277 

10.0111 

2.3280 

-1.1640 

3.4919 

2.5 

1.8273 

4.0945 

-5.9217 

2.26724 

10.0162 

2.3294 

-1.1647 

3.4940 

2.6 

1.8306 

4.0949 

-5.9255 

2.26425 

10.0203 

2.3303 

-1.1652 

3.4955 

2.7 

1.8331 

4.0947 

-5.9278 

2.26155 

10.0225 

2.3307 

-1.1653 

3.4960 

2.8 

1.8346 

4.0952 

-5.9298 

2.26055 

10.0249 

2.3311 

-1.1655 

3.4966 

2.9 

1.8365 

4.0949 

-5.9314 

2.25844 

10.0263 

2.3312 

-1.1656 

3.4968 

3 

1.8372 

4.0949 

-5.9321 

2.25775 

10.0270 

2.3314 

-1.1657 

3.4971 
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TABLE XVI: Simple 3D Hexagonal (Stacked Triangular) Spin Lattices with k = (1/3,1/3,0) and (1/3,1/3,1/2). 

Eigenvalues Aka and eigenvectors fi = fj, y , [i z \ in Cartesian coordinates of the magnetic dipole interaction tensor Gi(k) 
in Eq. (16c) for two values of the magnetic wavevector k in reciprocal lattice units (r.l.u.) and collinear magnetic moment 
alignments. The x axis is in the direction of the hexagonal a axis. The accuracy of the eigenvalues is estimated to be < ±0.001. 
Also shown are the differences between the eigenvalues for different ordering axes for a given k. 


c/a 

(U-o) 

[001] 

[100], [010] 

[001] - [100] 

(1 i I) 

[001] 

[100], [010] 

[001] - [100] 

0.5 

38.4703 

-19.2351 

57.7054 

-27.9927 

13.9964 

-41.9891 

0.6 

22.2841 

-11.1420 

33.4261 

-15.1774 

7.5887 

-22.7661 

0.7 

14.0899 

-7.0450 

21.1349 

-8.3821 

4.1910 

-12.5731 

0.8 

9.5488 

-4.7744 

14.3231 

-4.4347 

2.2173 

-6.6520 

0.9 

6.8755 

-3.4378 

10.3133 

-2.0181 

1.0090 

-3.0271 

1 

5.2388 

-2.6194 

7.8582 

-0.4912 

0.2456 

-0.7368 

1.1 

4.2100 

-2.1050 

6.3150 

0.4898 

-0.2449 

0.7346 

1.2 

3.5533 

-1.7767 

5.3300 

1.1260 

-0.5630 

1.6889 

1.3 

3.1295 

-1.5648 

4.6943 

1.5405 

-0.7702 

2.3107 

1.4 

2.8539 

-1.4270 

4.2809 

1.8130 

-0.9065 

2.7195 

1.5 

2.6742 

-1.3371 

4.0113 

1.9907 

-0.9954 

2.9861 

1.6 

2.5570 

-1.2785 

3.8356 

2.1078 

-1.0539 

3.1617 

1.7 

2.4793 

-1.2397 

3.7190 

2.1846 

-1.0923 

3.2770 

1.8 

2.4287 

-1.2144 

3.6431 

2.2347 

-1.1174 

3.3521 

1.9 

2.3953 

-1.1977 

3.5930 

2.2681 

-1.1340 

3.4021 

2 

2.3738 

-1.1869 

3.5607 

2.2898 

-1.1449 

3.4346 

2.1 

2.3597 

-1.1798 

3.5395 

2.3044 

-1.1522 

3.4566 

2.2 

2.3499 

-1.1749 

3.5248 

2.3139 

-1.1569 

3.4708 

2.3 

2.3435 

-1.1717 

3.5152 

2.3196 

-1.1598 

3.4794 

2.4 

2.3396 

-1.1698 

3.5095 

2.3244 

-1.1622 

3.4865 

2.5 

2.3371 

-1.1686 

3.5057 

2.3266 

-1.1633 

3.4899 

2.6 

2.3350 

-1.1675 

3.5025 

2.3283 

-1.1642 

3.4925 

2.7 

2.3341 

-1.1670 

3.5011 

2.3296 

-1.1648 

3.4945 

2.8 

2.3335 

-1.1667 

3.5002 

2.3304 

-1.1652 

3.4956 

2.9 

2.3329 

-1.1665 

3.4994 

2.3308 

-1.1654 

3.4962 

3 

2.3325 

-1.1662 

3.4987 

2.3311 

-1.1656 

3.4967 
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TABLE XVII: 3D Honeycomb Spin Lattices with k = (0,0,0) (ferromagnetic) and (1/2,0,0). Eigenvalues Aka and 
eigenvectors jl = [n x , fx y , fi z ] in Cartesian coordinates of the magnetic dipole interaction tensor Gj(k) in Eq. (16c) for two 
values of the magnetic wavevector k in reciprocal lattice units (r.l.u.) and collinear magnetic moment alignments. The x axis 
is in the direction of the hexagonal a axis. The accuracy of the eigenvalues is estimated to be < ±0.001. Also shown are the 
differences between the eigenvalues for different ordering axes for a given k. 


•^kct 

c/a 

(0,0,0) 

[001] 

[100], [010] 

[100] 

-[001] 

(i»,o) 

[#.M 

[ooi] 

[iXM 

-[001] 

[l 0 ] 

[2’ 2 ’T 

_ hd I ol 

1 2 5 2’ 

0.5 

17.8809 

-8.9405 

-26.8214 

-42.9938 

38.2651 

4.7287 

-33.5364 

47.7226 

0.6 

3.4722 

-1.7361 

-5.2083 

-31.1917 

21.8433 

9.3484 

-12.4948 

40.5401 

0.7 

-4.2063 

2.1032 

6.3095 

-24.5981 

13.3674 

11.2307 

-2.1367 

35.8289 

0.8 

-8.9442 

4.4721 

13.4163 

-20.5983 

8.5385 

12.0598 

3.5214 

32.6581 

0.9 

-12.1969 

6.0984 

18.2953 

-18.0445 

5.5981 

12.4464 

6.8483 

30.4909 

1 

-14.6077 

7.3038 

21.9115 

-16.3605 

3.7263 

12.6342 

8.9080 

28.9947 

1.1 

-16.4856 

8.2428 

24.7284 

-15.2276 

2.5001 

12.7275 

10.2274 

27.9552 

1.2 

-18.0090 

9.0045 

27.0134 

-14.4562 

1.6822 

12.7740 

11.0918 

27.2302 

1.3 

-19.2757 

9.6378 

28.9135 

-13.9270 

1.1274 

12.7996 

11.6723 

26.7266 

1.4 

-20.3528 

10.1764 

30.5292 

-13.5625 

0.7498 

12.8127 

12.0629 

26.3753 

1.5 

-21.2806 

10.6403 

31.9209 

-13.3097 

0.4905 

12.8193 

12.3288 

26.1290 

1.6 

-22.0891 

11.0445 

33.1336 

-13.1347 

0.3120 

12.8227 

12.5107 

25.9575 

1.7 

-22.8019 

11.4009 

34.2028 

-13.0134 

0.1889 

12.8245 

12.6356 

25.8379 

1.8 

-23.4352 

11.7176 

35.1528 

-12.9291 

0.1034 

12.8257 

12.7223 

25.7548 

1.9 

-24.0010 

12.0005 

36.0016 

-12.8708 

0.0446 

12.8262 

12.7816 

25.6971 

2 

-24.5114 

12.2557 

36.7671 

-12.8298 

0.0033 

12.8265 

12.8233 

25.6564 

2.1 

-24.9722 

12.4861 

37.4583 

-12.8016 

-0.0251 

12.8267 

12.8518 

25.6283 

2.2 

-25.3902 

12.6951 

38.0853 

-12.7818 

-0.0449 

12.8267 

12.8716 

25.6085 

2.3 

-25.7725 

12.8863 

38.6588 

-12.7680 

-0.0586 

12.8266 

12.8853 

25.5946 

2.4 

-26.1232 

13.0616 

39.1848 

-12.7585 

-0.0679 

12.8264 

12.8944 

25.5850 

2.5 

-26.4457 

13.2228 

39.6685 

-12.7519 

-0.0745 

12.8264 

12.9009 

25.5783 

2.6 

-26.7431 

13.3716 

40.1147 

-12.7475 

-0.0795 

12.8270 

12.9064 

25.5745 

2.7 

-27.0190 

13.5095 

40.5285 

-12.7443 

-0.0827 

12.8271 

12.9098 

25.5714 

2.8 

-27.2751 

13.6376 

40.9127 

-12.7421 

-0.0850 

12.8271 

12.9120 

25.5692 

2.9 

-27.5129 

13.7565 

41.2694 

-12.7402 

-0.0867 

12.8269 

12.9136 

25.5672 

3 

-27.7360 

13.8680 

41.6040 

-12.7392 

-0.0877 

12.8270 

12.9147 

25.5662 
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TABLE XVIII: 3D Honeycomb Spin Lattices with k = (0,0,0) and (0,0,1/2) (Neel-type in the ab plane and 
then also along the c axis, respectively. Eigenvalues Aka and eigenvectors fl = [ji x - ji y: // z ] in Cartesian coordinates of the 
magnetic dipole interaction tensor G;(k) in Eq. (16c) for two values of the magnetic wavevector k in reciprocal lattice units 
(r.l.u.) and collinear magnetic moment alignments. The x axis is in the direction of the hexagonal a axis. The accuracy of the 
eigenvalues is estimated to be < ±0.001. Also shown are the differences between the eigenvalues for different ordering axes for 
a given k. 


c/a 

(0,0,0) 

[001] 

[100], [010] 

[001] - [100] 

(0,0, i) 

[001] 

[100], [010] 

[001] - [100] 

0.5 

39.6849 

-19.8424 

59.5273 

-14.1112 

7.0556 

-21.1667 

0.6 

24.8305 

-12.4152 

37.2457 

-0.283961 

0.1420 

-0.425942 

0.7 

18.1343 

-9.0671 

27.2014 

6.17048 

-3.0852 

9.25572 

0.8 

14.9989 

-7.4994 

22.4983 

9.25024 

-4.6251 

13.8754 

0.9 

13.5052 

-6.7526 

20.2579 

10.7321 

-5.3660 

16.0981 

1 

12.7869 

-6.3934 

19.1803 

11.4468 

-5.7234 

17.1702 

1.1 

12.4402 

-6.2201 

18.6603 

11.7927 

-5.8963 

17.6890 

1.2 

12.2732 

-6.1366 

18.4099 

11.9585 

-5.9793 

17.9378 

1.3 

12.1920 

-6.0960 

18.2880 

12.0404 

-6.0202 

18.0606 

1.4 

12.1524 

-6.0762 

18.2286 

12.0792 

-6.0396 

18.1189 

1.5 

12.1341 

-6.0671 

18.2012 

12.0987 

-6.0493 

18.1480 

1.6 

12.1246 

-6.0623 

18.1868 

12.1075 

-6.0537 

18.1612 

1.7 

12.1206 

-6.0603 

18.1809 

12.1119 

-6.0560 

18.1679 

1.8 

12.1187 

-6.0593 

18.1780 

12.1146 

-6.0573 

18.1718 

1.9 

12.1164 

-6.0582 

18.1746 

12.1161 

-6.0580 

18.1741 

2 

12.1168 

-6.0584 

18.1752 

12.1160 

-6.0580 

18.1739 

2.1 

12.1162 

-6.0581 

18.1743 

12.1159 

-6.0579 

18.1738 

2.2 

12.1164 

-6.0582 

18.1746 

12.1166 

-6.0583 

18.1749 

2.3 

12.1161 

-6.0581 

18.1742 

12.1161 

-6.0581 

18.1742 

2.4 

12.1159 

-6.0579 

18.1738 

12.1164 

-6.0582 

18.1746 

2.5 

12.1162 

-6.0581 

18.1743 

12.1165 

-6.0583 

18.1748 

2.6 

12.1162 

-6.0581 

18.1743 

12.1167 

-6.0584 

18.1751 

2.7 

12.1165 

-6.0583 

18.1748 

12.1164 

-6.0582 

18.1745 

2.8 

12.1158 

-6.0579 

18.1737 

12.1156 

-6.0578 

18.1734 

2.9 

12.1164 

-6.0582 

18.1745 

12.1165 

-6.0582 

18.1747 

3 

12.1164 

-6.0582 

18.1746 

12.1164 

-6.0582 

18.1746 
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TABLE XIX: 3D Honeycomb Spin Lattices with k = 
(0,0,1/2) (FM intraplane and AFM interplane). Eigen¬ 
values Aka and eigenvectors (i = [fi x , fi y , fi z ] in Cartesian co¬ 
ordinates of the magnetic dipole interaction tensor Gi(k) in 
Eq. (16c) for the magnetic wavevector k in reciprocal lat¬ 
tice units (r.l.u.) and collinear magnetic moment alignments. 
The x axis is in the direction of the hexagonal a axis. The ac¬ 
curacy of the eigenvalues is estimated to be ~ ±0.001. Also 
shown are the differences between the eigenvalues for different 
ordering axes for a given k. 


c/a 

[001] 

^ka 

[100], [010] 

[100]-[001] 

0.5 

-47.097 

23.549 

70.646 

0.6 

-39.490 

19.745 

59.235 

0.7 

-36.499 

18.249 

54.748 

0.8 

-35.235 

17.618 

52.853 

0.9 

-34.674 

17.337 

52.011 

1 

-34.415 

17.208 

51.623 

1.1 

-34.295 

17.148 

51.443 

1.2 

-34.238 

17.119 

51.358 

1.3 

-34.211 

17.105 

51.316 

1.4 

-34.198 

17.099 

51.296 

1.5 

-34.191 

17.096 

51.287 

1.6 

-34.188 

17.094 

51.282 

1.7 

-34.186 

17.093 

51.279 

1.8 

-34.185 

17.093 

51.278 

1.9 

-34.185 

17.093 

51.278 

2 

-34.186 

17.093 

51.279 

2.1 

-34.185 

17.092 

51.277 

2.2 

-34.186 

17.093 

51.279 

2.3 

-34.186 

17.093 

51.279 

2.4 

-34.184 

17.092 

51.276 

2.5 

-34.185 

17.093 

51.278 

2.6 

-34.185 

17.093 

51.278 

2.7 

-34.186 

17.093 

51.280 

2.8 

-34.185 

17.092 

51.277 

2.9 

-34.182 

17.091 

51.274 

3 

-34.184 

17.092 

51.275 







